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ABSTRACT 

The analysis of scattered electromagnetic radiation is examined in 
terms of I'ie theory establishing the detailed form for the scattering 
functions. These parameters axe used to define the criterion i/hereby 
the composition and number densities of an atmospheric aerosol population 
may be assessed. Laser radar, or LIDAR, is the proposed remote sensing 
device and the governing system equations are developed. The problsm of 
data inversion is surveyed with emphasis on smoothing methods, statistical 
analyses, and iterative teclmiques. A discussion of the numerical 
stability of the solution is also presented. On the basis of the Mie 
model, inversion of the data provided by the LILilR probe is given as a 
rationale to effect an estimate of the true particle number distribution 
function in an atmosoheric cell. 
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I. riTRCDUGTIOI-y 


The advent of sophisticated optical devices in both the scientific 
and defense communities has pointedly demonstrated the need to better 
understand the primary propagating and interacting medium characteristics 
of the earth*s atmosphere. The intent of this paper is not to deal with 
the purpose of such devices, whether they are atmospheric pollutant mon¬ 
itors or infrared threat detectors. Rather the concern hers is the use 
of these devices in the real world environment and the limitations im¬ 
posed on them as a result of atmospheric optical signal degradation. 

This degradation is the result of manifold scattering effects such as 
Brillouin, Rayleigh, Mie, fluorescent and Raman scattering and includes 
the atmospheric constituents* differential absorption of optical 
radiation. CCl 

The atmospheric scattering and absorption centers include the air 
molecules, suspended particulates (aerosols) and water droplets in fog, 
rain or hail. The presence of aerosols strongly determines the trans¬ 
mission characteristics of the atmosphere and may even dominate the 
propagation characteristics# See Table I and figure 1 for a comparison 
of particle sizes and the range of number densities. 

At the small end of the size spectrum are the air molecules them¬ 
selves with radii of 10 ^ra and a concentration of 10*^^ per cm^. The 

large end includes rain drops with radii of 10 to lO^^m and a concen- 
-5 -2 3 

tration of 10 to 10 per cm , The **small** particles, defined by a 
radius much less than the v/avelength of radiation used as the probe, 
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Humber density of particles dll/dlog r (cm 


Average size distribution of nuclei in 
continental ciir masses 



Figure 1 
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obey the Rayleigh scattering model, exhibiting a cross section that 
decreases proportionally as the particle radius to wavelength ratio to 
the fourth power. It is this fourth power law for molecular scattering 
that Rayleigh used to explain the blue color of the daytime sky. 


TABLE I 


Particle t’/nes 

Particle Radius (/<a) 

Concentration (cm 

Aitken nuclei 

-3 -2 

10 ^ to 10 

10^ to 10^ 

haze particles 

10“^ to 1 

10 to 10^ 

fog droplets 

1 to 10 

10 to 10^ 

cloud droplets 

1 to 10 

10 to 300 


initial nuclei present in the atmosphere from the earth's surface 
[(J and/or meteoric sources. 

For ’.javelengthslarger than those of visible light, this scattering is 
very small. As the particle size approaches that of the radiation wave¬ 
length, the complete Mie theory is required to explain the scattering 
characteristics. The dependence on v/avelength is highly oscillatory and 
very complex. As the particle grows still larger, its behavior in the 
radiation field approximates that of a large solid object casting a 
shadow field in proportion to its geometrical cross-section. The Mie 
theory includes both Rayleigh scattering for small radii and the "non- 
selective" scattering in the limit of large particle radii. The theory 
is covered in some detail in the next section of this paper. 

Once one has both the atmospheric density and data on the aerosol 
composition and size distribution as functions of altitude, location, 
and time, the transmittance can be computed on a real-time basis by some 
rather sophisticated computer models. 

Aerosols are colloidal particles dispersed and suspended or falling 
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slowly in the gaseous mixture of the atmosphere* Some examples are 
smokes, haze, clouds, small droplet (less than 1/«m) fogs, and finely 
divided soils* These aerosols have as many sources as they have diff¬ 
erent forms and compositions* The chemical reaction of various atmo¬ 
spheric gasses whose minor components are nitrous oxides, terpenes, and 
hydrocarbons can produce solid particles* The presence of combustion 
products from inefficient burning, which may be locally critical, and, 
on a much larger scale, from forest fires and volcanic action serves 
only to increase the atmospheric load of particles that may take years 
to settle* Over large bodies of water, the dispersion of solutions into 
the near atmosphere region brings about the formation of fogs and rain 
as well as clouds* As a class, aerosols may either help or hinder the 
researcher* These particles may be used as tracers enhancing a received 
signal to detect atmospheric motions, distributions, and chemical com¬ 
position* However, they primarily provide a large, random source of 
environmental noise* In this latter mode, the particle’s presence tends 
to confuse the received signal at the sensor and provide noise that 
originates outside the system electronics, maJcing it difficult, if not 
impossible, to supress* ^6} Historically, aerosols have been divided 
into classes to • facilitate their theoretical modeling* These models 
involve several free parameters that are designed to ”best fit” the 
observed data* The classes are: 

1* Maritime - This group is characterized primarily of sea salt 
\vrith some water of hydration and a clilorine concentration slightly higher 
than sea water* Salt particle concentration decreases rapidly above a 
height of about 500 m* 
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2. Rural - The largest number of components are silicon, iron 
and some sulfates with approximately 10^ to 15 % organics. Of this total, 
about l/3 of the aerosol is hygroscopic, whereby air moisture is readily 
absorbed. This same class is also found at sea at higher altitudes vdth 
a laJTge particle density falling off rapidly as distance at sea 
increases. 

3. Urban - Combustion products and the issue from industrial 
processes maice up the bulk of this class. 

4. Troposphere - This is an atmosphere division lying above the 
surface boundary layer and consists of the iraral/maritime classes with¬ 
out large particles. 

5. Stratosphere - Primarily (30^ to 90 %) consisting of the 1-arge 
sulfate ions and particulates, this group also shows a periodic large 
influx during volcanic activity. 

The typical atmospheric analysis technique is to use two-ended sys¬ 
tems requiring a source of radiation of Icnov/n properties, a propagating 
path tlirough the atmosphere of specific length and reasonably loiovn 
conditions, and a stable receiver at the pathos end. Exp'eriments such 
as these are able to measure aerosol size and composition distributions 
to check out the various composition and scattering models. They also 
predict the performance characteristics of a particular optical device. 
However, two v;ell-separated platforms are required v;ith the attending 
problems of phasing and communications. Also, only low altitude, hori¬ 
zontal v;ork can be done. A one-ended system would eliminate the need 
for one of the platforms and enable a full altitude/azimuthal capability 
in a real-time environment. There is a system that can easily be 
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adapted to the one-ended process and is available today. The laser 
radar, LIDAR, system is similar in operation to Radar with optical 
wavelengths instead of radio waves. In fact, the acronym is somewhat 
the same, meaning "Light Detection and Ranging". Conceptually, a laser 
pulse is emitted into the medium to be analyzed. This pulse interacts 
v/ith the components of the medium and is both scattered and absorbed in 
distinct ways depending on the material type and density. Some portion 
of this pulse is returned as an echo to the point of origin where a 
telescope-lilce receiver is coupled to a photodetector that converts the 
received optical radiation into a proportional electrical signal. The 
pulse return timing is gated to the range of the atmospheric cell under¬ 
going analysis. The intensity of the echo, the echo's polarization 
change, and any measured frequency shifts all undergo involved processing. 
The outputs can consist of composition, type, and size distributions. 

Some state of the art applications of the technique include; metero- 
logical LIDAR investigations; middle atmospheric LIDAR studies; tropo¬ 
spheric chemistry and diffusion research; studies of atmospheric 
propagation and radiative transfer; absorption, Raman and fluorescent 
spectroscopic applications; and pressure-temperature measurements of 
atmospheric layers. A recent theoretical paper indicated that a 
series of measurements involving optical scattering and absorption can 
yield about 7 or 8 independent pieces of information concerning the 
aerosol size distribution within the analyzed cell, covering a range of 
particles with radii from 0.1^m to 5/<m. r2(3 

Consider such a single-ended LID;\R analysis system on board the type 
of mobile platform in use today by both scientific researchers and the 
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military. The optimum performance of a specific optical scanning/de¬ 
tection/tracking device could be evaluated in a real-time scenario vith- 
out the need to resort to generalized models applicable only in 
standardized areas. 
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II. GSIISRAL SCATTERING TIEORY 


Scattering is the process whereby an incident electromagnetic wave 
loses energy and some portion of this energy undergoes re-emission. In 
general, the scattered wavefront has a different propagation direction 
than the incident wave. See figure 2. 

For particles with a low number of electrons, the scattering can be 
computed using quantum mechanics. In this sense, the incident v;ave can 
only interact with the stationary states of the scattering particle's 
electronic system within a small bandt^idth of frequencies about the 
allowed transition levels. One can discuss Rayleigh scattering as the 
absorption of the correct energy photon and the subsequent re-emission 
of a photon of essentially the same energy. An energy level diagram 
would look like figure 3. 

Likewise, Raman scattering involves a definite frequency shift be¬ 
tween the incident and scattered radiation. If the scattered frequency 
is less than that of the incident, the process produces the Stokes 
lines spectroscopically. Aiid when the reverse is true, the anti-Stokes 
condition exists. C82 See figure 4* 

However, for molecular and larger particle scattering, the collection 
of electronic states becomes great enough to let us deal with the prob¬ 
lem classically. The particle is considered to be a solid dielectric 
with a complex index of refraction. 

One of the most widely studied problems in diffraction theory is 
electromagnetic and acoustic wave scattering by a sphere. The problem 
of diffraction of a plane-monochromatic wave by a homogeneous sphere 
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characterized by its radius and complex index of refraction vas rigor¬ 
ously solved by G. Mie in 190S, His analysis involved the solving of 
Maxwell's equations for the electromagnetic field when an incident plane 
wave interacts with a material interface. '.Then the discrete boundary 
separates media of different optical characteristics, a scattered wave 
is generated. 

Today, an extensive bibliography of approaches to this problem exists. 
The author will follow the analysis presented by Kerker (1969) 023 
throughout this section in order to present a consistent derivation. 

First, the salient features of the Maxwell theory will be listed, lead¬ 
ing to the wave equation. The solution satisfying the appropriate 
boundary conditions will be obtained; the resulting infinite series 
solution not only completely defines the scattered wave, but gives the 
electromagnetic states of the particules' interior as v;ell. The single 
particle scattering results 'will be compiled and extended to cover 
collective scattering. 

Max’well's equations, in rationalised MK3 units: MM 


Coulomb's law 

v-a 

= 

(1) 

jimp ere's law 

rxE 


(2) 

Faraday's lav; 

VxE 


(3) 

absence of free 

magnetic poles 

V ro 

(4) 

The sjTnbol denotes 

a vector quantity. 

Folloving are 

definitions 


of vectors used in derivation: 

- electric field intensity (volts/meter) 

^ - magnetic field intensity (amps/meter) 

2 

D - electric flux density (dielectric displacement, coulomb/metcr ) 
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^ - magnetic flux density (magnetic induction, weber/neter ) 

2 

J - electric current density (amps/meter ) 

- electric charge density (coulomb/meter'^) 
t - tLne (seconds) 

These equations will define the electronagnetic field config^J.ration 
for all points in space and uniquely detemine the field vectors given 
the media characteristics belou, 

a = eg (5) 

^ ( 6 ) 

I- ^ ^ (7) 

The media parameters thus defined are: 

iy - permittivity (electric jrductive capacity, farads/meter) 

- permeability (magnetic inductive capacity, henrys/meter) 

^ - conductivity (specific conductance, l/ohn-meter) 

In free space (vacuo): 

3.8542 X 10“^^ farad/m 
4'’n{ 10“’^ henry/m 

As a consequence of the invariance of the Maxwell equations, the quantity 


^6)”~ has the units of velocity and, again, in free space: 


C= X 10' m/sec 


6 


(S) 


This last quantity is the speed of light in vacuum, ’./hen the quantities 
ando'are independent of direction, the region in space they 
characterize is said to be isotropic. Two core equations ere required 
to complete the picture. 

Lorents Force £ = X §) (9) 

I'e’,rton's second law of motion ^ (13) 
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'■rnere: ^ - force on charged particle (newtons) 

- velocity of charged particle (ui/s) 

^ - nass of charged particle (kg) 

2 

a - acceleration of charged particle (m/s ) 

/>/ 

All macroscopic, linear, electromagnetic phenomenon are described by 
equations (1) tlirough (10) vrlth the associated quantity definitions. 

The classical d^^-Tisrics of interacting charged particles and electro¬ 
magnetic fields are thus specified. 

A certain degree of convenience is obtained using the definitions of 
the tvro following auxiliary vectors: 



(n) 


(12) 


Tnere: ^ - electric polarization (c/m ) 

- magnetic polarization (V^) 

One can see that, in free space, these quantities vanish identically and 
they may then be considered as a measure of the effect of natter on the 
local fields. Some dimensionless quantities may also be defined: 

- §/e, (13) 

'.•"cere; \(^ - specific 6 (relative dielectric constant) 

- s?ecificy^(relative magnetic permeability) 

The media may further be described with: 

05) 

06) 

I'/here: electric susceptability 

magnetic susceptability 

The previous set of equations can be combined to obtain: 
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Xe = Ke -1 07) 

Xm= Xm-l (18) 

Maxuell' s equations apply to regions v/here the aedia may be repre¬ 
sented by the continuous functions and O', ’.fnen a medium discon¬ 

tinuity is encountered, such as the conditions across a regional 
interface, the fields must satisfy a specific set of relations at the 
boundary. See figure 5. Let; 

2 

^5 - surface current density on interface S (A/^^ ) 

- surface charge density on interface S {c/a) 

The Maxuell equations (l) - (4) are cast in point form for each 
differential point in space. One can apply various integral theorems 
to change the same equations to integral form and use them to define the 
conditions that the fields must satisfy by crossing the interface.0'^ 
The boundary conditions obtained are inherent in the Maxwell equations. 

Boundary conditions: 

(a) the tangential component of ^ is continuous 

(£2 ~ ^ £ =0 (19) 


(b) The tangential component of K undergoes a discontinuity equal to 

a= Is (20) 

(c) The normal component of D, undergoes a discontinuity equal to 

2 = ^5 (21) 

(d) the normal component of ^ is continuous 

(g2-g.)-R. = O (22) 

For all cases of dielectric scattering, it is assumed that there is no 
free surface charge or current density. This assumption begins to brealc 
down for large particles in the presence of storm activity. However, ;;e 
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interface 3 


Figure 5 
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’.■all, for simplicity, assume that ^ = 0 and(^^ = 0, Therefore, the 
tangential components of E_ and H and the normal components of D and ^ 
are continuous across the boundary. 

Charge conservation is included in the Maxwell theory and is written 
as: 

continuity equation (23) 

Vihere: - volume current density (Vm ) 

- volume charge density (c/m ) 

Throughout the follo’wing derivation, a few vector identities will be 
required, 

ITote: continuously differentiable vector field 

0 - continuously differentiable scalar field 

V • (V X A ) = o 

V = vV 
r x(^i 2 f) = o 

= A • {V0) + ^2^(v-A) 

V X (vxA.') = vCv-d) - v^/i 

If the fields described by (l) - (4) are invariant under a gauge trans¬ 
formation, then they are v;ell defined by the preceding formalism. 

VJe define: 

vector potential field 
0 - scalar potential field 
in the following relations. 

B - VXA ( 25 ) 

Compare this to equations (4) and (24). Using equation (25) in (3), we 
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get the following: 

^ (l. = O 

How, refering to ( 24 ), 


=> E +• 




■i= -vjzS' 


, thus 


E - 




■4-V^ (26) 

A and 0 are not unique. Let g be any scalar function, then the new 
potential can be defined as: 

A' = ^ - VS (27) 

0'r j2)' + ^9 (28) 

Now the fields E and 3, are recomputed, 

0'=rx (d-^3) = B- Vx(r9')'B 

=> § 

£'=‘itC4-va)-v(^+-|r3)= '^A-'!'sf+ir(’’3)-r 




The physical observables are the fields E and ^ and it can be seen 
that the gauge transformation defined by (27) and (23) leaves these fields 
unchanged. Thus, the Maxwell equations are gauge invariant, 

V/e now turn from formalism to a more practical aspect of electro¬ 
magnetic theory and introduce the Wave Equation: 

r X ( 6 ): VX (r X E) i ^ (V X B) 

Note that the time and space operations of differentiation are commutative 
for these continuously differentiable functions, 

(0, C 2 w): r (v-£ ) - r^£ = '-|r 

The assumption is always made that all space and time derivations of the 
parameters 6,^, anddare either slowly varying \d.th respect to the rate 
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of change of the fields themselves oc identically equal to zero. 


=> V rv-e)- 


a), (2), (5), C^) : V. <:££)= £V-£=/’ 

=> V • E = /Vt 



The same argument also applies to the quantity/^/8 concerning its rate 
of change with the rate of change of the fields. 



(29) 


This equation, (29), is the uell-kno’.jn nonhomo gen eous equation for damped 
wave notion - the '.lave Equation. SLmilarly, V X (2); 


rx (vxfci) = 


V • ^ \i - O 1=> V • H = O 



(30) 


The relations for non-dissipative media are obtained by setting o' - 0, 
Both equations (29) and (30) actually represent 3 equations each. Each 
component of the vector field must satisfy these equations for the vector 
field itself to do so. Let the symbol u correspond to any of the six 
components of the two fields E, and The scalar wave equation can be 
expressed by: 



(31) 
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The fields described by (31) represent a transfer of energy in both space 
and tine. The electromagnetic wave has a time periodicity that sug¬ 
gests a fora for the energy flux crossing a unit area per unit time: 

S = E H (32) 

Equation (32) is the Poynting vector equation, whose tine averaged 
magnitude over the period of the wave yields the wave intensity, 

I = '/r/jll dt (33) 

whereTis the period of the wave. The field symbols E and H in (32) are 
the real parts of the complex fields v/hose conponents satisfy (31). 

The form of the scalar wave equation suggests a harmonic time de¬ 
pendence of the form: 

U (iC.'t) = Hz) (34) 

\jhere o) is the angular frequency, i =V(-1), and r is the usual three- 
dimensional position vector. Note that: 




LujU 


_ Z . I 

and = u) LL 


d-t 

then, (31): V% - c co U ^ co ^ U = O 

=> V e cr u = O 

Define the propagation constant k by: 

j<2 ~ l/A CTud 


(35) 

then: V U = O 

Equation (36) is the Helmholtz Wave Equation, teother form of (35) is: 


K = — i, /3| 




(37) 

(38) 

(39) 


It is the/3, term that determines the damping undergone by the wave per 
unit path length; the energy is dissipated as Joule heat in the medium 


26 

























v/henCS' >0, 


Recalling that the quantity 


V - 


(40) 


is a velocity which, in free space is c, the relation for an electro¬ 
magnetic wave in a nediun described by light speed V is 

^ (41) 

where: X - wavelength 
■0 - frequency 
V - velocity of the wave 

Also, 2Tr\) = u) (42) 

and this is the angular frequency, u> , mentioned earlier. In free space, 

(43) 

The index of refraction, n, is defined as 

(-U) 


=c 


n = c/v 
,): n 

=> n X = Aq 


(41), (43), (44): n= = ^°/x 


(45) 


It is clear that A.© is the wavelength when the v;ave is propagated in 
free space. In the sane nediun. O' = 0 implies that 

/3, =1 0 => '-''I = to - u>/c 


=> Ko= ZV-A 


(46) 


This is the free space propagation constant. 

The complex refractive index, n, can be defined as 


m 


= k/Ko= n (I - iK) 


(47) 


where: hC - index of absorption (extinction) 

As an example, suppose that there is a plane v;ave propagating along the 
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z - ajcis whose x - component 

= /I* < 


(47): 


Luoi 


- A& & 

=.>E^-Ae 

=>£;/ - Ae 


iuot -inf^e> O ~ i X) z 


The aotenuation of the v;ave is seen e:’:plicitly as e 


-rf' X: 


■in ejroansion 


of this enample shows that the resultant plane wave solution leads 
directly to the Stokes paraneters. Assusie the propagation direction is 
s, as stated earlier, and that the nonoclironatic wave is traveling 
unbound in a source free isotronic nediun. It is assumed that the wave 


by completely specified in 

the ::y-piane 

v:ave. Thus: 



Aik . AJk 

h X ^ y 


(34), (36) 


- O 

ioit 

S =0 




- ae 




Substituting a trial solution, 

=> ak^e‘’'‘+ K'-ae = 

=•> f(z) = 

From the definition of k, -icz is selected as the vaJ-ue. 

^ - '■ 

=:y U 12:,t; = ae 

This equation describes the fora for each of the six vector components 


of E 

and H. 

Recalling (l 



d£<i 


(1) 

dx ^ 


^ 3z. 

But, 

ax = 

o) “ 

0 '> 


£ 

1 

3 


= o 

'dB.z, 


and U £(x ^ ^ 

£7, = o 

B i(^u)i-/<z) 

= a y e 
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( 3 ), ( 6 ) 


vx(r= = —y , 


- 'aT-^=|7^ > -0 => 

It is quite sufficient to describe the behavior cf^ as the value of ^ 
can be readily obtained using equations (l) through (4). Given the form 
of u(s,t): 

3 H 

a-t 

IL 

52 : 

=> 

=> i u> 

.;oT:e: =■ <-tj/K 

The term/^^/k can be reduced using (35) to the folloid.ng: 


"> 6t 



3 

My 
at ’ 

Z i cO 

II 

1 

3 

-~y 

a z ^ 

= -ik 

-t i k. 


11 

X 

^ 1 

- 



J< _ 

^ CO ^ 


^(t - 1 

-> ^ 

/^UO 


L <J 

K 

^ UJ 

^^'to y 


Defining the intrinsic impedance of the medium for plane naves as 




y? 


(4S) 




-E 


-V 




Hy = 

linally, 1.-°- plane wave is: 

t. (a)t - Kz. ) 


r Ae 

= Be‘ 

£,= 0 


i Cct)± - Kz. +c^-■) 


(49) 


where and ofiare arbitrary phase angles to account for the general his¬ 


tory of the wave. The !I v/ave is also determined by equation (49): 

t _ i(u)t-Kz foTj) 

^ X “ ~To ^ ® 

LCcOi -K2, 




(50) 


= O 

;/e see that: 


E-H 
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(51) 


=> E • H = O 

Thus, T and "" are *oernendicular to each other for the olane wave 

solution. 

The concept of polarized light can novj be introduced uith the help 
of equation (49),^",, and^. The general plane wave solution describes 
an elliptical rotation of the T vector in the olane pemendiculax to the 
direction of propagation. Defining the phase difference as«/: 


^ ^ Sy - S, 


(52) 


The representation of the state of a polarized wave can then be uniquely 
deternined using the Stokes parameters, which are: 

= fi^ t sf (53) 

s, = = So cos 2 v* cos 2 X ( 54 ) 

2 f\b C.05 if = Sq sm z'f'Cos z X ( 55 ) 

S3= 2A6 SINcTc So sin Z X (56) 

The lengths and 3, and the angles andX are illustrated in figure 6. 

The quantity (cot - kz) can be eliminated from (49) using equation 
(52), where E and refer to the real parts of the complex v/ave: 

This is the equation of an ellipse. The ellipse in figure 6 is drawn 
"right-handed", that is, the head of the electric field vector foliov;3 
the arrows when viewing the wavefront face on. In order to exhibit right- 
handed polarization, sinS>0 and 0^X4fl"/A* Linear polarization exists 
when: 

S = j= 0,11, i 2, t 3, - • (58) 

and circular polarization corresponds to 

fej-r, j=tl,±3,ts, ±7, ... ( 59 ) 
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Pictorial definition of the 
quantities ^ 



2A 


ta,. X = I V«, X « %), {04 TT) 


Pigure 6 
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The Stokes paraneters satisfy: 

^0 = ^2 + ^3 ( 60 ) 

and a complete state of polarisation is obtained from only three in¬ 
dependent parameters. 

Natural light is defined as the state where the intensity of light 
is the same in any direction perpendicular to the propagation vector, l£. 
The phases of the various components in natural light are randomly 
varying so that no correlation exists between them. 


=r> 


n « 


= S 2 = S3 = O 


I - - 7 - J ^ — W C 6 I ) 

The combination of eUiptically polarized beans can be accomplished using; 


So = .^ Sq\ j S, - ^ S,i. j S2- ^ ^3 ■ ^ ^31 

Then for non-identical polarization states for all beans: 

sj > s? ^ ^ s" 


( 62 ) 


(63) 


Partially polarized light can be described as a superposition of natural 
light and a set of incoherent beams: 


and 


sr, , ^ 3 - S 3 , 


Finally, the fraction of polarization, F, can be defined as the intensity 
ratio of the polarized part of the wave to the total wave superposition: 

^ ' ^o/(So+So) (o^F^l) (65) 

It can be shown that the parameters of a scattered wave are linearly 
dependent on the parameters of the incident wave. This relation 


is written as follows: 

- 


13 = (i/R‘;5i-, 


( 66 ) 


where: R - observation distance to scatterer 
i - incident wave 
s - scattered wave 


32 

















I = 



The nodified Stokes paraneters are defined first by writing equation 


(49) ! 

is: 

il^it 

LU)i 




Ex 


= £2® 


(67) 

then, 

I. - ie, 




(63) 






(69) 


U = I R 




(70) 


V 3 2 r 

'”1^1 ^z*} 



(71) 

The 4X4 Stokes matrix o, looks like 

this: 




ip„r 

IfJ' 


■Lm 


C7 = 

lU' 




i(72) 


2 Reffiifai) 





2 Ir»nCf|(f2|) 

1 : 


Re 



v;here the f^^. are defined in terms of the scattering functions given in 
1^2- For spherical particles, the quantities fj2 = ^21 “ ^ ^ 

contains four constants. 

Having laid the foundation for the Kie treatment of scattering 
phenomenon, let us now consider a plane electromagnetic wave incident on 
an isotropic, homogeneous sphere of arbitrary size. The object interferes 
^d.th the free propagation of the wave and secondary iraves are generated 
that are predictable by Maxwell's equations and the boundary conditions; 
and these secondary v;aves form a scattered radiation field. The incident 
wave should be linearly polarized for two reasons. First, if the light 
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is unpolarised, then it can be resolved into tuo linearly polarised 
conponents each of which act independently of the other. One only needs 
to consider the interaction of the sphere with a plane (linear) polarised 
wave and then apply the principle of superposition to obtain the total 
effect. The LIDAR laser pulse is plane polarised, by definition, and 
this is the wave-type that is considered for tropospheric probing. A 
geometric diagram of the problem is shown in figure 7, suggesting a 
spherical coordinate system as the most convenient for analysis. 

The fields of electromagnetic energy in the region of space filled 
by the particle and the incident wave are resolved into three conponents: 

a. the incident wave - E., H. 

b. the oarticle interior wave - E , H 

'■r^ >rr 

c. the scattered wave - E , H 

For these quantities to properly define the fields, they must each 
satisfy Manwell’s equations (l) - (4), the IJave equations (29), (30), 
and the boundary conditions (19) - (22). The internal field E , H 
must natch the external field S. + E , H. + H in accordance uith 
these relations. 

KerkerQ^et al., suggests that the derivation deal vdth the scalar 
wave equation, (31). Two new functions may be introduced, the electric 
Hertz vector, and the magnetic Hertz vector, Ti^, and are defined as 


follows: 



(73) 

(74) 

(75) 

(76) 
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Scattering of 


a plane wave by a sphere 


r7 
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incident plane v/ave 


Fig’ire 7 
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Suitable transformations may be applied to define other Hertz vectors 


as long as the fields are left invariant. The vectors thus defined must 
satisfy the vector wave equation: 



(77) 

(78) 


where the^ and result from distributions of electric and magnetic 


dipoles. Therefore, two vectors can be used to describe all the fields 
present. 

The transverse magnetic ’.;ave (TM), also called the electric wave, 
has a zero value for the radial component of its magnetic field intensity, 

- 0, Likewise, the transverse electric wave (TE), the magnetic wave, 
has a zero radial component of electric field intensity, - 0, The 
Ti! wave is conceived as being the result of oscillating electric di¬ 
poles in a spatial region, and the TH wave as the result of oscillating 
magnetic dipoles. The total solution is a construct of these tv;o 
wavefronts. 

The Kertz-Debye potentials are: 



(79) 

(80) 


These potentials are the solutions of the scalar wave equation. Solving 
(79) and (SO) for the potentials and applying the appropriate boundary 
conditions, one can then find the field vectors describing the TM and 
TS i;aves. The total field may be derived by the direct addition of the 
component waves. 

The follo’wing equations for the total field vector components can be 
obtedned by applying these concepts to the above relations. In spherical 
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coordinates, we have: 

Er -- ^ + K^CrT?:) + 0 

^e~^l 6 *^ 2 e =” "arT© 1^2 rsin© 7 ^ (''"z') 




Efi = E ■*■ rsin e “ ^ 2 ,"r © C'" 

Hr - H,r +l-l2r " ^ ('■‘n'j)+ K^CriTi) 

TTsTo© (r /7^) + r ^r9© 


H;:^rHlgf + H2/= K,-j:-|^(rTr>Vi:4;7r© 

vdth the propagation constant defined in the follo\7ing fore by (35): 

1:^ = " K, K’2(S7), V7hereK|=- LoJC' + o' (SS), K 2 * ^ dJ (c9) 

yt< has been dropped as all nedia are assuned to be non-nagnetic. 

Recalling the scalar wave ecuation for a function vrith harnonic 


(31) 

(32) 
(83) 

(34) 

(35) 

(36) 


u. j , iwt 

iine depenoence, e , 

iujt 


we nave: 


( 36 ) + K^a = o 


(34) U---P(r)£“ 


The spatial part f then satisfies: 

2 


K -PCc) - O (90) 

The Ilertz-Debye potentials are constracted such that they are solutions 
to tliis ;:ave ecnation. Let IT (r, 0 ) “ f(r) represent oither 

potential, 7r, ov\-. 

V^TT -h K^iT o (91) 

This equation is the one that must be solved; in spherical coordinates: 

1 f 1 ^ c)1T\ J- </_ z . . 

r C^T*’ ^ *0 r*sir»© de ( d O K 'TT = O 

It ca:i be solved using the technique of separation of variables, whore 
fTis assuned to be a product function: 

Tr(c) = RC^) ©Ce) (93) 


-> C'"^) n O 5 r R 


(q/\ 
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=> ^(^'“* 11 ) = («; 

=> ® '°°' 
(92! beooMs;^ ^ I 5 (siM S If.) 4.-^^ . -R^R®! 


' ir.GG ' T7 ' ":, then -jo r . ultr ; 

_r ^(‘rR . J i_AiN6l@’)., 


r 


1T : 

J_ S I. ^ 2 

~h~^ - “K r 

i^+KV^ P_A/swe^J_l^n 

— "*’|®^'N© V d©y'*5siM^S ^ 0^_j 


R 3r^ ® SIN © ^©V 

r 


(97) 


rRCr) = O 


R ari 

Tl'_Is equation is equal to a. constant which nay be cheched by tal:ing the 
derivatives of both sides with respect to r, © , . Tliis constajit has 

the forn of n(n + I): 

r R + = nfn + O 

Thus the radial equation is: 

d 

Using (99), equation (97) becones: 

_j_ e h®) j _ ,h 

®SIN©^( reJ'^^siM^e = -r\(n + i) 

LW^ Jq{ n(n^i)s.M 0J= ^ 

By the sane differentiation aTgu-ment, each side of this equation is 

2 

equal to a constant, m : 


uaT^ion is: 


(98) 


Thus, the 0 -angle eouation is: _ 


(99) 


Using (98) and (99), (97) talces the forn of: 

1 ^ 

Thus, the 0-angular equation is: 


( 100 ) 


In the above equations, n nay talce on any integer value and n nay be 
any one of the follovring set: -n, ..., 0, ,.., n. Equations (99) and 
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(100) conbir.e to forn the Spherical Hamonics vrhich.are defined: 


n < 

/here A is a nornalization constant to ensure orthonornality of the 


n 


function set of Y. 

The solution to (lOO) is a straightforuar’d coabination of siiies and 

2 . 


cosines as a is al’.'ays positive: 

o.„ C05 + Ism sm Cm 


( 102 ) 


lcu 3 .tion (99) is the faniliar .'.ssociative Legendre equa.tion having 
the ren-L-an solut: 

(103) 


®(e) = 


(rn) 

wliere the (cos © ) ore the associative legendre ool^rorniohs. 


Jsiiig Rodri^aes^ fornula : 

q ^ 


(1C4) 


where :: = cos ^ and the negative values of m are included, -n^ra ^ n. 
■ihen n - 0, then (99) is the legendre equation whose solutions exe the 
lefqsndre nol'^aionials: 

'<d\^ 


Lcuation (IO4) may be rexaritten as: 


(105) 

( 106 ) 


The complete spherical harmonic solution includes the normalisation 
constant and reuTites the azimuthal solution as a comple:: e:n)onenticl. 
These functions are completely orthonormal as well. Thus: 

Y^Ce, Z (-0'"/^ pt\co. e)a (,07) 

The (-1)*^ term is a phase factor defined as the Condon-Ghortly phase, and 
may, in some texts, be dropped. 

(hote: in the remainder of this paper, I './ill use ^ instead of zeta. 
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uMch is used in the literature.) 

The Radial equation, (9£), has solutions defined in terns of the 
Bessel functions: 


(lOS) 

XnCX'') = "yf Kr A/n^'/2.CKr) ( 109 ) 

=> r = Cn CxO X< (Kr) (no) 


These are the Pdcatti-Bessel functions, t-zhere J 1 =. 
integral order Bessel and heumann functions. By using 
definitions of the Bessel functions, one can ’..rite the 
function in the follou-ing derivative forn: 

v/here 3 = x + iy 


II . X are the half 
n+v 

the infinite series 
Ricatti-Bessel 

(m) 

( 112 ) 

(113) 


The function: 

-- fn X„CKr) = .JFy (,, 4 ) 

vanishes in the linit as kr approaches infinity, and v;ill be useful in 

(2) 

later work. The function H\/(kr) is the half integral order Hankel 

n+T 

function of the second kind. 

How it can be seen that the solution to the scalar wave equation is 
characterised by the nunbers n, m: 

rJrr’= rR„(r)©'"\e) 

and ( 102 ), ( 103 ), and (llO): 

=> rfr^"''’(r,e»= (115) 

The general solution nay now be obtained as a linear superposition 
of these particular solutions. 

rTTCO =11 rlT^^^r.e.C*) ('I6) 

n:o mi-n 


40 









Equation (116) is the general solution of the scalar wave equation in 
spherical coordinates* 

The fields in and about the particle can be described in terns of 
their individual pairs of potential functions; 
a* the incident wave - Tr,‘, n 

b. the interior wave - Tir, 

5 5 

c. the scattered wave - rr^,rr; 

The sphere, whose origin coincides with that of the spherical coordinate 
system and is of radius a, is isotropic and homogeneous, optically 
described by a complex propagation constant, : 

(47) ~ ^0 (117) 

The isotropic, homogeneous medium is a dielectric described by a real 
constant, k 2 * 

(113) 


The relative refractive index is defined; 

K,/k^ ^ no, /noi (119) 

The incident plane wave propagates along the positive z-axis and is con¬ 
sidered to be of unit amplitude with polarization ajcis parallel to the 


x-axis* 




( 120 ) 


In order to facilitate matching the boundary conditions for the Debye 
potentials on the surface r = a, (120) is written in the form of (Il6), 
realizing that the function Xn(^2^^ becomes infinite at the origin 
where the field (l20) does not: 




( 121 ) 

( 122 ) 
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The unperturbed incident wave is now well defined. Because the incident 
wave propagates along a coordinate axis (a), parameter m = 1, As a 
result, the Debye potentials all have the associated legendre poly¬ 
nomials of degree 1. The internal particle fields must be of the same 
form, xd.thout ^^(k,r), and must also be finite at the origin. Also, the 
arbitrary coefficients must exist in order to properly map the conditions 
at the boundary; 

C.„'''r,CK ,0 P„‘'’(C0S e) tos fS (123) 

rTT/ r Tf “i i"" dn'f’aCx,'') s) s.« (124) 

The scattered wave must vanish at infinity without entailing an infinite 
energy and the function ^r,(kr) contributes just this property while 
maintaining the proper form for71^(r). Arbitrary coefficients are a- 


gain necessary: 

^3 _L 

= -7- 


K’-z- -C 


n-i 


, n(n!o Pn''^ (coi e) COS 0 (l 25) 

'■’i'N-ih I-/'"' 

The boundary conditions at the interface ensure the continuity of the 
tangential comnonents of E and H. In soherical coordinates, this trans- 


lates into: 

Eq (r = a) = Eq C r = a ) 
E(zf c rzCL) =• £ ( 2 r Cr= a) 
He ( r-a) = ( r-a) 

(r=o.) = H^Vr=a) 


(127) 

(123) 

(129) 

(130) 


where (l) refers to the sphere interior and (2) the exterior. The mix¬ 
ture of terms in TTj and at first glance malces it difficult to apply 
the boundary conditions. However, by forming the appropriate linear 
combinations of (31) - (36), the boundary conditions can be recast into 
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an equivalent fora for the Debye potentials. Define: 


^ (siuseg) * 4 ? 

(131) 


(132) 

o'-is-Hs (siNsHii^) 

(133) 


(134) 

Then the boundary conditions (127) - (130) become: 

A^^'^Cr^a) = A^^\r-ia) 

(135) 

B"’(r = o.> = B'*’cr = 0 

(136) 

C^'^Cr-.a.)= c‘^\r:a) 

(137) 


(138) 

'lote that: — TT 

(139) 


(140) 


The quantities A and C renove r "iTj, and 3 and D remove rW,. The decoupled 
equations then result in boundary conditions on the Debye potentials. 


Substituting (31) - (36) into (l3l) - (134) yields: 

rA (Ul) 

rB = COS (U2) 

/y^i ( 143 ) 


For equation (135) - (138) to be valid for all 0,on the coordinate 


surface, r - a, requires the various bracketed arguments in (I4l) - 
( 144 ) to be equivalent between regions (I) and (2). Thus: 



= [lcrrr.)y;i 

(145) 

[^r 

- [i?crn;)ty 

(146) 

[K r<r,tL 


(147) 


= r>C.rrr(]-\ 

(143) 


43 













■Jou, using (139) and ( 14 O): 



(149) 


(150) 

YCi'Vir,'' = ir,^) 

(151) 

= 'l<'jV(lThv3 

(152) 


Equations (149) - (152) are the required boundary conditions on the 
Debye potentials. Substitute (121), (123), and (125) into (I5l): 


K 


<0 *• 


(0 “» 


Since all terns in the series expansion are linearly independent of each 
other for each value of n and for all &, p(; then the corresponding terns 


in the series nust be equal: 

0 ) . . sx-'W 


^ '^n ^ [r„ (K^ci)- c,„J„ 

Cr,y„ («,<=-)= (itih - a n ^ n (K I'^O 


K* yC] 

The quaiitity^ can be reduced by using equations (87), (89), and, 
Ki 

recalling that the frequency of an electromagnetic wave is unchanged in 


traversing the interface: 

- k!’ k7 kT X? , 

k’ k,‘" X T K':’ ' Xz’ ' ■ 

Applying similar argunrents to the remaining boundary condition yields a 

Ki m.Ko / 

term of the form: -- = —J-rr- = m,/rn, =: rn 

tn,,Ko ' 

The four boundary conditions then become: 

Cn Y'n (153) 

’inY/Cx.'^)- 'rnCr„'(HzCi)-):.n^n(KiO.)] ( 154 ) 

Cr,YnCK,0= (155) 

= C'^„(KzO-yhn'^n(K^CL^ (156) 

Bote; Y'(^r)=-^y(.Kr) 

This set of equations uniquely determines the coefficients a^, b^^, 
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However, only those corresponding to the scattered wave are 


and d , 
n 


inportant here. Define the following paraneters: 

=: k^a -- = zrr 


ti - K.a, = ivKjCl - 2’rr^,^ - 




(157) 

( 158 ) 


where is the vacuun wavelength. 

The quantity «x is called the diaensionless size paraneter. Solve 
for c from (155) and use the result in (153). Also solve for d from 


n 


(156) and use the value in (154). He obtain: 

% = (•< > Yn'C^i -/« /"/I 

b _ 

' ^3^ c-<') V'A -'fn (/3}3/; C’^) 

All the field vectors are now uniquely described given the Icnovna 


(159) 

( 160 ) 


parameters ando<. This data enables the coefficients a , b , c , and 
^ n’ n’ n’ 

d to be computed and their potentials determined, Finally, the vectors 
n 

are obtained from (Si) - (36). The formal solution is complete. 
Practical light scattering observations are performed in the far field, 
or wave zone, Tliis condition is such that where n is the 

Ricatti-Bessel function’s order parameter. In tliis approximation, the 
previous relations become simpler. The as^miptotic e;rjan 3 ion of the 
Hanlcel function becomes []} : 

( 162 ) 

The fields becone transverse due to the radial conponent decay as 
(A/r) compared to the (A/r) dependence for the non-radial components. 
Finally, using (Si) - (£6 ) vdth (125), (126) and (l6l), (l62 ) for the 
scattered wavs ortly ( the tends 'to zero for large 'c^r with re¬ 

spect to ^^(h^*) a-s does ) yields: 
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E0 -- 


mj. 

-H0’ 


- -S. 


& ’‘f<,r 


s/V v5( C©) 


v;here the annlitude functions nre: 

an+/ 


(163) 

(164) 

(165) 

( 166 ) 


S, C®)= 7;^ [a„7r„('a$9)tLnX,Ccose)] 

and tho an^'ular functions are: 

ir„ c©) = 

r^ffcos e) = ^ P3’(«se) 

The energy flow in the scattered wave can be calculated by using 
Fomenting* 3 theoren: 

S = E^ziHs) (169) 

where - represents the complex conjugate of the quantity. The intensity 


(167) 

(16S) 


[“unction is defined as follows: 


: (e) = |sr©)| 


( 170 ) 


The intensity of the scattered radiation for unit incident intensity 
polarised in the ^ and© direction are: 

1© " (171) 

5^^ z / • 

l0 = q Sl^ pc, (172) 

2 

Note that the intensity decays as l/r as it should for this spherically 
spreading wave. The scattering plane is defined as that plane which 
contains the incident direction and the direction of the scattered wave 
(©,< 20 . Then i^ is the wave perpendicular to the plane and i^ is the 
wave parallel to it. The phase relation betv/een E 0 and S^being arbi¬ 
trary is the condition that defines an elliptically polarised wave. In¬ 
tuitively, one can think of each polarised component coming from that 
diroction as being inherent in the incident wave. Therefore, the 1015 
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caused by an incident uave of intensity cos ^nclarized parsdlel to 
the scattering plane. 

The phase difference that exists betv;een these two parts of the 


scattering wave is: 

tan Re(Sz) +ImCS,)I,virS 2 ) 


(173) 


IJnoolarized incident light gives rise to a scattered bean of intensity: 


= SirV - (‘•'^^0 

with degree of polarisation given by: 

i. I - L 2 


P = 


Cl + '•2 


(174) 


(175) 


Fixing the scattering plane in space to be the ys-plane, all ob¬ 
servations by definition tal:e place in the scattering plane. The in¬ 
cident electric vector has its direction in the xy-plane at an angle X 
to the y-azis. The x-a::is is designated vertical and the yz-plane 


horizontal. This is shorn in figure 8, Thus, the vertical ('/) and 
horizontal (H) components of the intensity in the scattered wave for 
unit incidence become: 



^ XvCX) 

-it 

II 

pi 

^ X 



(176) 

I® 

= rH(5c) 

- JL 





(177) 

The Stoke 

s parameters 

; are: 






So 

- ■ c 

- L/rP^r^ V 

i, 5 /A/ 

X -h c 

; CCS 

^x) 


(173) 

S, 

- r; 

, 3/A/ 

" X - C2 

CCS 

^x) 


(179) 

Sz 

“ y h’-i. 

3///' X 

CCS 

X cos 

cf 

(130) 

Si 

- / 

^ Z 

3//V X 

C03 

X S//V 

/ 

(1S1) 


’Then discussing cross sections, one is speaJcing of the amount of 
area that the particle effectively presents to the incident beam by its 
ability to spread out the incident energy. The complex index of 
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Scattering plane geometry 



Figure 8 









refraction includes absorption o.s well as scattering. Define the ex¬ 
tinction cross section as the sua of the cross sections for scattering 
and absorption : 

(1S2) 


By integrating the real part of the tiae averaged Poynting vector for 
the total external field: 

I = ^ (1S3) 

the total outward flow of energy is obtained. For unit incident energy, 
the various energy losses represented by the integral resulting from 


above are the desired cross sections: 


^SCQ, = -^ ^^1 ] 

cre^Lt 


(154) 

(155) 


The efficiency factor is defined as the cross section per unit 
2 

geometrical area (ffa ) of the particle: 

Q= (186) 

Terms of the form^ X becomesince«< =ETrOL/x, Thus: 

Ztt' CO 

Qsco, = ^ ^ rant/') (187) 

Qext = ,;:fr+ [Re (138) 

Ilote that '^sco. and are independent of the state of polarisation of 

the incident wave. 

Kerker et al., gives an excellent review of the difficulties 

involved in the numerical computation of the previously listed scattering 
functions. The functions needed to be evaluated are a , b , Tf , . 

Some discussion of the equations is given in the Appendix. 

Using approximations for small arguments, the- fourth power law 
attributed to Rayleigh scattering can be directly obtained. The first 
three scattering coefficients are: 
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0.2 


-- ^/3 C 

/ mJ: --L) 

L z / 

; 

/5 <■ 

/nililA o< 

C i/h’* +3 J 

“ MS i" 



(1S9) 

( 190 ) 

(191) 


plus terms in higher powers For the Rayleigh criterion«x<<l, 


only a^ 

contributes. 

From 

the Appendix: 



rr, ( 

005 9) 

- 1 



(192) 


X/ 

(105 ©) 

It 

o 

S 0 


(193) 

(165), 

( 166 ): 

s,(e) 

_ 3 

~ t 

--•s- 

a, ;i, 

(cos ©) 




s.(e)^ 

- 3 

z 

d, r, 

Ccos © ) 


(170): 


i,(^) 

_ 9 

~ V 






(i-C®) 

_ 9 
“ V 


cos^& 



For unpolarized incident light, for example, sunlight; 

(174): lu =S^ ^ la.,l"‘(H- cos^©) 

_ cjn-w j I m^-l ' 

(189); la = SzrP^ T 

)\\ 

1^(1 + co5^ ©) (194) 

This equation accurately predicts the blue sky phenomenon as the shorter 
wavelengths scatter much more than the other visible wavelengths in 
daylight. Equation (175) yields the observed degree of polarization: 

P-. 

1 Id, 1^ + )(t,r cos© 1 

p. I - cos^ 0 (195) 

I +C05 o 

A plot of this function on a polar graph indicates a maximum of 1, 
complete polarization, at a scattering angle of 90°. This has long been 
observed in the day sky. The error associated with the trioncation after 
n = 1 is 5^ if a/X is kept 0,05. This limit is the required Rayleigh 


(157): Xu 


A X 


_ 817-'^'' 


c><- ( I + C(3S ’• ©j 


criterion. 
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The cross sections and efficiencies represented by equations (184.)- 
(l8S) give the total integrated differential cross section with the 
integral formed over the entire solid angle sphere about the particle. 
The differential cross section is: 


d cr ( a) 
d-0. 




( 196 ) 


where: ^ - number of photons scattered into a solid angle per unit 

time per particle 

^ - number of photons crossing a unit area per unit time 
dn - Si-fi 0 d& df/ ^ the differential solid angle 
The intensity is obtained by multiplying the number of photons at fre¬ 
quency "j) by Planck's constant, h, as hv> is the energy per photon. 

h = b.lD2(?2 X 10 Toul<i/s ^6 (197) 

As written, equation (196) describes a scattered energy per incident 
intensity. The scattered energy becomes a scattered intensity when 
divided by the distance to the observation point, squared, for then the 
solid angle becomes the subtended area, 

d A ^ r ^ d IT- (198) 


utilizing the above, 

dl ^ 


dA 


v/here dl represents the scattered intensity over dA per 
intensity. The equations describing such an intensity, 
may be cast in the form: 

H fr 

where 0 (d) = C05 ^ or Sin ^ 

i ©) = or C ^ 


(199) 
unit incident 
(171) and (172), 


( 200 ) 
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I is the dl per unit intensity 


Thus: 




( 201 ) 


This equation, when integrated over all angles &, ^ yields the cross 

section given in equations (1S4) and (l85). Hereafter, the synbol \/±ll 

aean dCJ/d-aas defined by (201), This is in keeping with the notation in 

the literature* 2 . / n / \ 

=;^ i ((/) 


( 202 ) 


This area represents the scattered energy per particle per unit intensity. 
The theory aay now be extended to multiple particle systems vdth the 
addition of the following definitions. Let n be the total number of 
particles in a region of differentiable volume, cross-sectional area and 
extent given by dV, dA, and dx respectively. 


Then: 

dV = dxdA 

( 203 ) 

and 

n = NdV = HdxdA 

( 204 ) 

where H - particle number density 



Assume the scattering to be incoherent among the particles. 

Then 

the total O' for all the particles is 

just 


= n O' 


(205) 


where ^ - scattered energy per unit intensity for one particle 
The amount of intensity scattered into dA is the scattered energy 
divided by dA. 

O' n _ scattered intensity (206) 

“3^ " unit incident intensity 


.'mother useful parameter is the scattering coefficient: 


/3 


scattering intensity _ 

incident intensity-path length 


(207) 
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’.rhere tlie nath length d:c is that traversed by the incident radiation. 
Ccnbining the alcove equations: 




=> (5 = (yN 


(203) 


Until nov7, r represented the radius of the radiation field and a the 
particle radius. Again, to be consistent vith the literature, r must be 
assigned to the particle radius. Using the new notation, if there exists 
a population of scattering centers each described by U(r), then the a- 

?rom all the ^articles of radius r between r and 


r 4* dr is: 


£ 7 ( (3 — cr (^ni ^ ©, 0^^ c(A/ Cr) 


(209) 


The quantity dU(r) may be WTitten as: 


( 210 ) 


dU 


'.■7here ^ is the particle number density per particle radius. Then: 


<^0 = o- (^ ) dr 


( 211 ) 


( 212 ) 


Let the lower and upper radii limits be r^ and r 2 respectively; vre 
have: , ,, 

In the case of o' then/3 -> the extinction coefficient 

for a volume of scattering particles defined by dV and equation (212) 
would represent the total energy loss by an incident beam passing 
through a volume of depth dx. Assuming an x (path length) dependence 
for , then the amount of intensity lost per unit incident inten- 

drex) 


= "/3(xjdx 

■X 


rex) 

XCx) ^ e)cp 


(213) 

(214) 
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where I^is the incident intensity to the volume at x-0. Tliis last e- 
quation is loiovm as Bouguer's law. [ 23 ) It describes the degradation of 
the incident bean as it traverses a material path. In order to express 
the scattered intensity into a solid angle dii^we recall that the in¬ 
tensity scattered per particle into solid angle dH was given by equation 
( 200 ), Let the parameter r in that equation be replaced by R as r now 
describes the particle radius. 

(215) 

As before, the number of particles in the cell, II, applied to equation 
(215) yields the total scattered intensity per unit radii per volume; 
then, for all particles of radius r between r^ and 

^ ^ I (e) g C^) dr 

=> I = r ( nn, e, 

lor unpolarized radiation: 

c (no, e) 3 (0) = ^ 

where i^ and i^ are given by equations (l65) - (I 68 ) and (l70). For the 
polarized radiation: 

i (m,oc (-2 (218) 

Equations (212) and (2l6) are the specific results of the theory 
used to investigate aerosol parameters and v/ill be discussed in the next 
two sections of this paper. First, some comments on the collective 
scattering phenomenon are in order. It is v;ell knovni that the more 
ordered the scattering centers become, the more coherent the scattering 
due to the phase relations between sites. m This leads to a reduced 
baclcscatter ;a.th corresponding increase in forward scattering. In 
crystalline substances, there is essentially no backscatter, when off- 


( 216 ) 


(217) 
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resonance; and the incident beam merely undergoes propagation velocity 
changes. The condition met by a gas at atmospheric pressure or less, 
where the number density is low on a radiation v;avelength scale, ensures 
that the scattering is incoherent, allowing the summing of intensities 
from all the scattering centers. The Mie theory for a collection of 
particles requires a random distribution of scattering sites to achieve 
this incoherent intensity accumulation. Quantitatively, the near-field 
interactions between the particle may be ingored if the average distance 
between sites is at least ten times their radius. This criterion is 
satisfied even for very dense fogs and the independent scattering model 
computations are valid for practical tropospheric scattering. However, 
significant number densities ^ri_Ll lead to a large flux of previously 
scattered light in the region and add multiple scattering effects to the 
received signal, A quick checl: on the validity of the single scatterer 
model is to verify linearity between the observed effects, I and/3, and 
the number density of scattering particles. If the relation is non¬ 
linear, one must suspect the existence of coherent and/or multiple 
scattering. Alternatively, if the relationship is linear, it is certain 
that the independent, single scattering conditions are present, A 
measure of this condition may be obtained from equation (214), (2^ Let; 



(219) 


be the optical death of the medium, then: 

-r X 


X C^) = XqS 


( 220 ) 


The conditions are: 

rx<o.i (transmission > 9055 ) - single scattering predominates 
o.i^ Tx < 0.3 (transmission 75 to 95^) - multiple scattering 


effects may be removed via correction to the 
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single scatter model 


T'y (transmission <15%) - the complete, non-ana- 
lytic, approximate multiple scattering theory 
applies 

Virtually all aerosols in the troposphere satisfy the first criterion. 

The literature abounds with material containing graphs of the Mie 
scattering functions for many specific cases. The Appendix outlines one 
such case for a specified set of conditions, Kowe'v^er, some general 
information can be reduced from this plethora of data. The amount of 
forv/ard scattering is about equal to that of the backscattering at a 
size parameter i^) equal to 0,01, Asciis increased, the amount of for- 
v;ard scattering also increases with a liJce decrease in the backscatter, 

1, the forward scatter is on the order of 100 times that of the 
backscatter. Increasing still further, there e:d.sts a unique scattering 
maximum with the existence of smaller subsequent maxima and minima. 

The scattering function increases smoothly as (Rayleigh) up untiled 
^0,5* Thereafter, the structure is highly dependent on the complex 
index of refraction. The large oscillations in the behavior as®< is 
vaTied become smoothed as the absorption index increases. As the num¬ 
ber density, as a function of the particle radius, taices on mors of a 
detailed structure, there is a further smoothing in the scattering 
functions as small differences tend to be averaged out. Conceptually, 
one has at hand either an actuaJ. particle distribution or a specified 
model of such. 

Given IT(r), for each r, we can compute the required I^('x',©) and sum 

j. 

the total intensities using 'l(r) as a weighting distribution. This 
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total l(e^, 0) then appears as the snoothed f’onction, 

Liice the scattering function publications, there are also a large 
quantity of analytical nunber distribution nodeis in use. The verifi¬ 
cation and application of those nodeis has been the primary purpose of 
□uch of the aerosol e:q;erinents conducted in the past 20 yeans, and the 
nor:: is continuing today. Is conputers beccne faster, real tine inver¬ 
sion analysis ‘./ill lead to a nore accurate determination of the "true'- 
nunber distributions, Tvo such nodels shall be nentioned here, 3one 
natural aerosols nay be described with a n;;o-paraneter function utiliairn 
c* V0 ccci.3^"b j Cj 3 '''d c s'*’ s.”"r'^ ss 


presented by Junge Q*] ; 


cl ^ 


( 221 ) 


d'h'r^ . 

■ I I ~ < c! o o 


previously defined in (210) end nay varp' from a.bout 


2,2 to 4.0* As an enanple, for silicates in the 0,1^n to 10,0range, 
one nay use; 

r, = O.OM/um , y V =n3 

Another such model is in realitp'’ a set of nodels that includes the Junge 
model as a special case, Deimendjian [l] has proposed: 


fCr) = 


c^^/ (f) ^ -br 

~TZ — - 0.r e. 


( 222 ) 


(note; the ^ used here is different from the size parameter,) The 
versatility of this approach is the presence of four adjustable parameters 
(a,<9<', b, y) sJ-lov/ing for a vide range of experinental data to be fitted. 
Three specific cases of this model are; 


Haze C (continental): -P T 


Haze II (maritime) 
Cloud 


- b'Tr" 


■P re 

o 4 -1.5 r 
^ r®e 
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Host of the other proposed nodal distributions utilise soae forra of 


equation (222) as their basis £und canned routines for their conputer 
modelling are readily available. 

T'jo other effects that have been omitted in the present theory are 
Doppler shifts and Brillouin scattering. The presence of random 
(thermah) anoustic vaves in an atmospheric cell will add a degree of 
order whereby Brillouin scattering of the incident '.rave can occur off 
of the resultant density fluctuations. However, this effect is extremely 
•;ecL; for gaseous materials and the degree of coherence added to the 
structure of the gas is generally undetectable. Similarly, this motion 
of the gas will cause Doppler broadening of the wavelength of the in¬ 
cident radiation adding to the natural line v/idth of the source, Tven 
for the sharpest sources, la-sers, the natural linewidths are on the 
order of several ilHs, ’.rnercas, for almost all atmospheric conditions, 
the e:q:ected Doppler shifts sne on the order of a few tenths of a Mlis. 
Therefore, the sliift is generally lost in the natural linevridth and the 
relative motion between the detector and the atmospheric cell may be 
neglected. 


III. LID.4R 

The sophistication of the laser systems available today, coupled 
with the fast computer processing that can cheaply be obtained, has 
singled out the laser radar, LIDAR, system for increased use as a tropo¬ 
spheric e:<perimental device. The LIDAR should have the ability to per¬ 
form a wavelength scan with selectable bandv/idth receivers. It is the 
atmospheric attenuation of the laser pulse and the introduction of 
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random phase and amplitude changes that combine to form the scattered 
pulse, or echo. The subsequent inversion of this data then allows the 
e:q3erimenter to assess the optical characteristics of the atmosphere 
on a real-time basis. With the system geared to perform azimuthal and 
elevation scanning as well as range gating, one can then map out the 
atmospheric optical behaviour within the visible hemisphere. The two 
most important atmospheric parameters to be considered for a pulsed 
LIDAR system are the extinction and scattering coefficients. The total 
extinction coefficient may be uritten; 


(3 - /3^ ^ /3 A t 


(223) 


where R - Rayleigh (molecular) scattering 
M - Hie (aerosol + other) scattering 


A - Molecular and Aerosol absorption 
0 - Other processes (Raman, fluorescence, etc) 


The other processes tahen together are typically two or more orders of 
magnitude wealcer than the Rayleigh effect and tend to get lost in the 
signal noise, 

Rayleigh scattering is obtained from the Hie solution assuming «>c <<1 
using equations (184) and (1S9) - (191): 



(224) 


Rayleigh scattering for air molecules assumes m is real and m^:^1; 


^ = ZZa./x: 



(225) 


where V = (4/3)fra^, the volume of one 


the volume of one molecule. Let N - l/V, the 

s 


molecular number density, then: 



( 226 ) 
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where u = (6+3^)/(6-7(f), the depolarization factor, required (227) 
to accoimt for the anisotropy of air nolecules, for which cT = 0,035. 
"he quantity 

afe d'i - r\ nun'll /"ij-rs Q<ii> nr^e. 

(228) 


^(e,g^) = 0.0b07l (1 + 0.932 cos ^9) 


also accounts for the anisotropy. Let (°cu be the air nolecule number 
density (vice II), then: 


^ 229 ) 

where it is assumed all air molecules have approximately the same radius 
of 10 ^^m. 

The Mie extinction coefficient is dependent on particle radius 
through equation (212) whose only unknoi-m is dll(r)/dr. The unique sharp 
v;avelength feature of the LIDAR allows one to work well away from the 
characteristic absorption lines of the sample. Conversely, it also 
allows one to type the scatterer composition. For example; 

Laser 1: A.,®^ ^ absorption (resonance) line of a particular 


species; 


—> /v3| = /3p^i + (^/r>i ^ (^jqj 


( 230 ) 


where is the resonant absorption extinction coefficient wiiich is 


strongly dependent on X ’with tabulated and 

= <^A 

_ f 

The is the number density of the species. 
Laser 2: A.^jwst off the resonance; 


(231) 




flZ 


Due to the relatively wealc X dependence for 

/S, - /3^ a /s;' 

For an analysis cell a distance R away: 


( 232 ) 


(233) 


6o 




(2u) I(r) = 

and bade at the receiver: « 

I'CO.KR'le 


Then; 


x; 


lo, e'^-^ 


•oZ 


Thus, the f®* for the particular species nay, 
this equation using the inversion techniques 
For the typical nagnitudes of the considered 


(234) 

(235) 

in principle, be found fron 
outlined in the next section, 
cross sections, see Table II. 


TABLE II 


Process 

o' (cm /ster) 

Mie 

10 -2’ 

10“^ 

Rayleigh 

10-27 


Raman 

0 

-P 

0 

1 

0 

10“^ 

Resonance Raman 

10 “^^ 


Fluorescence 

4 10“’^ 



Sinilar to radar in operation, the LIDAR has a perfomance equation 
v:hereby the system characteristics and the appropriate atmospheric para¬ 
meters vhcih contribute to the signal are related.[ 23 ] Consider a 
pulsed LIDAR operating in monostaticniqd.e (receiver/transmitter co-located) 
being used to analyze a specific atmospheric cell via its optical prop¬ 
erties. The problem geometry is outlined in figure 9. In this figure: 

hj. - cell height 

R - slant range to cell 

n -receiver solid aagle 
r 

^ - transmitter solid angle 

u 

L - cell length ’’gated" via receiver on-time 
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Monostatic laser radar 


/ 



laser 

transmitter 


laser 

receiver 


4- 


(exaggerated) 


Figure 9 
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A - cell area seen by receiver 
r 

- cell area seen by transmitter 

The xollo'/ing derivation shall consist of counting photons at each 
stage along the beam path; converting to intensities has been previousl 
discussed. 

1 . I'lx - number of photons transmitted per pulse 
Vj. - photon frequency (transmitted) 
hv>^ - energy’’ per ohoton 
"xhV, - energy per pulse 

- laser outout oouer 


- oulse width 




(236) 


2. for Goncstatic LIDAH, the scattering angle is IcO degrees (ff} 
for backs ca,tter. Assume that all oackscattered oho tons are received: 


=> 




(237) 


The cell volume is determined by tlie Lim.'H timing sequence and is: 

(23s) 


The solid angle; 


n = A/R' 


V = LAta* 


( 239 ; 

(240) 


The pulse width constrains the ma:-cimum system resolution by limiting the 
minimum cell length achievable: 

(241) 


•“min = 


where c is the speed of light in the medium. 

This condition assumes tha.t the pulse tail refloctixig off the front face 
cell does not interfere 'Tith the pulse front reflecting off the 
=1 1 
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3 


. I;. - nunber oT photons incident on the cell 

(Sx - atmospheric extinction coefficient for the given 


transnitted 

N)u =^(6 ■-C%h)d5 


(242) 


Ho'/ever, to account for beam forming requirements, assume an optical 


system efficiency of 




(243) 


4. The baclcscattering coefficient nay be defined in accordance 
with equation (212): 

^ a ( n )~ ( 3 ( 1 T ) =, r,© = / (244) 

v/here (’a’(n) is the usual notation. This equation represents the back- 
scattered intensity per unit incident intensity per unit solid angle per 
unit path length. The signal returned from the cell is due to this 
process alone, Bouguers’ lav (214) states, for PaCn) assumed constant 


throughout cell length: 


X-e 


^ (ocr (J7)L 


(245) 


’•/here I is that oorticn of the incident beam that nasses tlirough L un¬ 


affected, The scattered mortion is I = I - I : 

^ s 0 




(246) 


where IT - the nunber of uhotons backscattered and a ra.tio of intensities 
s 

is the sajae as the ratio of the number of photons. T^gical tropospheric 
aerosols yield: 


fO crCn)L « I 

(if = 0.02, overall error 1(4) 


(247) 


The restriction this daces on 


1 j -Up chi 


r. irrcer limit of several hundrec 


meters; Thus: 

1.1 
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-/°<T(rr)L 


- e Ci. (0<f (ir) L 

:> n' = A/l P(r(fr)L 


(24G) 

(249) 


9>.e 'orinie on indicates the number of ohotons ocr solid anrie that have 


been oaclcscattered. Combining equations (243) end (249): 

^ C^)L e%p I'f [3^ els \ 


( 250 ) 


5 . low tho echo must pass through the sane propagation path on its 
return to the receiver as it did on its nay to the cell. However, the 
extinction coefficient nay be different, due to the accumulation of 
23 roce3Ees lil:e Raman scattering, frillcuin scattering, Doppler siiift, 
etc., though individually negligible. 

- number of photons at receiver mirror 

(239), (250): 


(251) 


j (250): r T1 , 1 

= Nt a-(n-)Lfl. etc L WtM 




( 252 ) 


6 , The optical system vri.ll collect the echo photons received at an 
efficiency of ^ and convert them to an electrical signal via a detector 
v;ith a quantum efficiency o£^ : 


f = Nnn 


(253) 


where li^ - number of photons received and counted per transmitted pulse, 
due to bachscatter from the atmospheric cell of interest. 


Tnus: 


Mr = ^cr rrr)Lft.-^r7q^e,Xpi~I d^] 

I - I ■ ■ " ■ ' ^ ^ 


(254) 


This equation nay be recast in the form of intensities if the photon 

population returned is large enough that individual photons need not be 

counted. Let D be the optical system aperaturc (D, - transmitter; D - 

"0 z* 

receiver), then the intensity becomes: 
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(255) 


I 5 NK v>/d 

Re^^riting equation (254): r r'Vu./x /o . D 1 

T It ^ ^^2lds3 

Ir - -—--- -4 


( 256 ) 


Thus, the important parameters for the LIDAR analysis of an atmospheric 
cell are the scatterer number densities and associated cross sections, 
lor a given system design, one has available the quantities it>%- 
A^, D^, D^, and can select the values for , L, and R, If one 
assumes a slow moving target cell, then the pulse round trip time, t^, is 


tp = 2R/c 

R.= Cc./0^r 


(257) 

( 258 ) 


Thus, one determines the range, R, by turning the receiver optics on a 
time, t , after the oulse is emitted. In order to avoid the oroblem of 
range ambiguity, lilce radar, one must ensure that t^ be less than the 

X 

period between pulses: 

tr < l/PRF (250) 

v/here PR? - pulse repetition frequency 
The unambiguous range limit may be written; 

= (c/e)/PRF (260) 

where R^ - unambiguous range 

The cell length, L, is determined by the time that the receiver is turned 
off again a time At after t^. For the entire cell to be scanned, the 
pulse must traverse L tv;ice in At; 

L= C c/2) At (261) 

and the receiver must be turned off at time; 

t - tp t A t (262) 

One can then select L and R by programming a timing sequence for At and 
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t^. Finally, the interpulse period, T, which is a system deadtime be¬ 
tween pulses available for processing, becomes: 

T. (I/PRF) - r* (263) 

Nov; one measures and '))f and applies equation (256) to determine the 
desired quantities ^’it) andThis inversion is covered in detail 
in the next section. Note that -when the equation for CJfrom 
(202) is used to replace the ®,in (244), then equation (256) 

for has the correct form as given by ( 216 ) for the scattered 
intensities. 

Thus, a model has been constructed ’whereby the system parameters are 
related to the cell unlmovwns in a definite manner. But first, some con¬ 
sideration must be given to system noise. During the time the receiver 
gate is open, t, the system vri.ll collect and count the background sky 
photons as vjell as the signal. The continuous ’’dark current” in the 
detector yields an electrical noise as a result of dark counts. These 
counts arise from photoelectrons being emitted from the detector photo¬ 
cathode surface even 'without an illumination present. Each of these 
noise sources may be reduced by proper design of the equipment and ex¬ 
periment. In the case of a detector, its dark current, I^, is either 
specified or may be measured in the laboratory: 

^ ( 264 ) 

v;here - number of dark counts per sec 

e - electron charge per coimt 
G - detector count gain 

The number of dark counts detected as part of the return signal in a 
time At is Nj^;4t, and : 


r r-j 

o / 




Mj) - (265) 

Let N„ represent the background count rate such that lUAt is the nun- 

D " D 

ber of background counts received per laser pulse* Assume a bacicground 
spectral brightness, that accounts for all the natural and artificial 
bacicground sources. This V/-^ represents a power density per solid angle 
in the background sky. Then: 

^ ^ ( 266 ) 

h V 

where ~ previously defined 

V ~ frequency of background photons (atA) 

A A - receiver banduith, centered about 'X 
Let 11^ be the total number of noise counts per laser pulse; using 


equation (26l) for At: 

Mr, = (2L/c)(Mj,tNs) (267) 

Due to the pulse repitition frequency, FRF, the noise count rate: 

CPRF)Nn (2^> 

v/here N - noise count rate 

The statistical nature of noise yields an effective RIIS noise count rate 

of'/if. The signal has a count rate of: 

Sr = (PRF)Nr (269) 

v;here S - count rate from signal echo alone 
r 

Thus, the total counting rate is: 

S = PRF (Nr N) (270) 

Talcing the RII3 of the signal plus noise, one can form the signal-to- 
noise ratio: 

siH = 5r/^^3- (271) 

Assuming a one second time averaging: 





































(272) 


S/N = V PRF 

If the noise average, :I , is not acciirately neasured, the II. is replaced 
by 211^ (measured) in the above equa.tion. Ilote that: 

1 . 0 => (s/A/)o ^ ^CPRF)Nr (273) 

2. if Nn= - C^l^)(i/\rr (274) 

Therefore, the S/lI PR? and by increasing PRP, one increases the signal- 
to-noise ratio, Kov;ever, via equation (260), the unambiguous range de¬ 
creases faster than the increase in S/lI. Since a lange S/li is desired 
uith good raaige resolution, one selects a detector v/ith a lou j^7^(-'Q) and 
designs the system for lov; bachground reception (P^). If the time aver- 
aging had been performed over a period of T seconds such that "RT T 


pulses had been transmitted and received: 

S/ri = Vprft i^r/v/Tw/r 

lonbining equations (265)-(267) and (275): 


N = 


sTpRF r A/r 




(275) 


(276) 


Inamination of the above equa.tion and (254) leads one to some general 
system noise qucl-ifications. 

1. One desires a lo'u detector I^ and a high''^^. 

2. Select a narrow receiver spectral bandwidth,A X; another reason 
for the choice of a laser. 

3. Reducing L reduces tho noiso as w'ell, but the transmitter must 
put out greater powci- to maintain the sax’.e signal. Then, for tho same 
average po'wer, the PR? must be lowered and the effect in S/lI is only 

\|PR? as compered to the 1 decrease and lij. incroase, increasing S/’’', 

In general, high po’wer/low PR? yields a better S/lI than a low power/ 
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4. A snailrequires a lov;er bean divergence at the transnitter 
and on optical matching 'oetx-ieen the receiver and the echo bean pattern 
(another laser capabilit;/). 

5. Without any noise, there still exists a finite (S/:I)q vliich 
represents the "shot noise" statistical nature of the signal. It must be 
emphasised that all the quantities in equation (276) represent the nesn 


os tne .inancc.'oet s'as’saoj-os. 


6 . Mr « Mr =i> S/N ~ Nr 

If Mr Nn => S/^J ^ V777 

The actual tropospheric neasurenent will include both 
(noleculan) and hie (aerosol) scattering. As an exanple: 
a. single ’.rater droplet, r = 1 (hie) 


:277) 


(27S) 


lO”® cm 


(^ '^ 10 cm 

-> (^m 


~3 


t a”7 _ “I 


10 


cm 


"" 10 


*2& E 

cm 


e 

=> 0R 


\'l - 3 

10 Cm 




Thus, the addition (3(^ to the total(3^ is on the order of a far; per¬ 
cent and any quantitative nea.s’ur9nent must t?J;e into a.ccount both the 
aerosol and noleculan scattering. The t’virbidity is defined as the ratio 
of the total scattering Intensity to the Rayleigh sca-tterlng intensity: 

= It/Tr Xzb) 

oone researchers attempt the neasurenent of this quaaitity independent of 
tile intensity, then use the res’ult to deduce the aerosol I... .Another 
nethod assumes that the a.erosol concentration above a certain height 
(e.g. ^ 30 I:n) is sufficiently lov; to allou it to be neglected. Then 
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(220) 


equation (256): 


[Text, one computes from (256) assuming the conditions of the U. S. 

Standard Atmosphere and equation (226). ThuLs is the e:cpected echo 

intensity due to molecular scattering alone. The ratio betvreen the mea- 

sured and the expected yields a scale factor, b, that corrects the 

standard atmosphere for the local conditions. Using this factor, one 

assumes a density distribution for the local atmosphere the same as that 

for the standard and commutes the emected value of i_ at the cell's 

R 

location. This intensity is then scaled using b. The measured value of 
I_ is then corrected for this Ip(b) and the result becomes I,., the Mie 


intensity: 


(2S1) 


For other means of separating the molecular and aerosol intensities from 
a measurement, the reader is referred to the ;;or!-: by Hall, et. al, Lvl 
This Rayleigh-atmosphere assumption is the starting point for the in¬ 
version process. 

The Standard Atmosphere may be represented via the folloving equations 

temperature T('’F) = 59.31 - 0.00357 (2S2) 

pressure P(a4m)= i.022. ezp(-^.157/ (2S3) 

Mass density _ O,002MZ5ey:p (-3-3^5^ (284) 

(slugs/ft-^) 

where h may vary fix)n 0 to 37,000 feet.[)3j The more precise form for 
these relations are also represented as canned computer routines. 

The critical term in equation (256) is the v/here the sub¬ 

script M refers to the aerosol distribution alone. The extinction com¬ 
ponent in the exponent of e in this equation represents the scattering 
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and absorption cross sections integrated over the entire 4'ff steradian 
sphere. If one assumes some representative model for the aerosol dis¬ 
tribution, as was done for the Rayleigh case, then the only unloiovm in 
(256) is The two averaging processes, integrating aver 4 'n' and 

integrating along the path R, reduces the sensitivity of the exponential 
term in (256) on the actual aerosol distribution. After completing the 
inversion process, one may input the best estimate for tliis distribution 
into the exponential term in (256) and then repeat the inversion. Ko’-p- 
ever, the processing time involved may prove to be the final constraint. 
The LIDAR experiment yields the scattered Hie intensity from the 
aerosol population in the analyzed atmospheric cell as a function of 
the pre-selected parpmeters at the laser (e.g, wavelength and polarization 
scanning). With this data as input, the system is prepared to perform 
the inversion routine. Some computational needs are reduced for the case 
of laser backscatter as the scattering angle remains fixed at ISO degrees. 
This baclcscattering condition simplifies several calculations as out¬ 
lined in the Appendix, 


IV. MATHEMATICAL DR/ERSION 

Tlie assumption made throughout this paper implies that the e:<peri- 
menter need only apply the exact Mie theory to his observations and per¬ 
form a simple data inversion to discern the aerosol properties of 
interest. However, this inversion is by no means simple. Aerosol prop¬ 
erties may be separated into two distinct catagories. The first con¬ 
siders those parameters that are explicit in the previously defined 
scattering functions; i.e,, spherical aerosols, aerosol sizes, low 
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number densities, and the relative indices of refraction. Then there 
are those aerosol characteristics that require some approximations (in 
many cases, severe ones) to the theory; i.e,, non-spherical shapes, 
anisotropy of the shapes with respect to the scattering plane, and high 
number densities leading to mioltiple scattering effects. The first few 
LIDilR experiments investigated the actual aerosol ranges for these para¬ 
meters, They found that in excess of 90% of the real aerosols studied, 
the spherically symmetric single scattering theory was valid under favor¬ 
able atmospheric conditions. These conditions require the absence of 
dense water clouds in the atmospheric cell and a depolarization of the 
incident beam limited to a ma:d.mum of 2%, with ^,5%> of less actually 
measured.[^ 3 ] In the following analysis, these favorable conditions 
will be assumed, as will the requirements of single particle scattering. 
These conditions are enhanced by the common field of view shared by the 
LIDAR transmitter and receiver. 

Recall equation (216) for the scattered intensity a distance R av/ay 
from a cell containing an aerosol population described by dll/dr with 
radii limits r, to r^: - / ui, 

,=.« 1 144. " 

The first simplification is to replace the particle radii continuum with 
a set of distinct particle classes each with a characteristic radius, 
number density and complex index of refraction. Thus: 

= *“■> n > ''iv ') •>•••) A, (235) 


v/here m. - j 


th 


oarticle refractive index 


.th 


r^. - j particle radius 

N. - particle number density 
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X - incident radiation v/avelength 
^ - scattering angles 

The large number of unknovnis that need to be determined will ensure a 
non-unique inverted solution from the limited sample of observations-.a- 
vailable. There is also the added difficulty of variance in the sup¬ 
posedly Icnovai data v/hose uncertainties at best generate unbounded noise 
in the final results. Therefore, some attempt must be made at reducing 
the initially large number of unlmao’.,Tis. 

The dependence of the intensity on the size parameter, has been 
discussed previously. Briefly, the non-monotonic fluctuations as is 
varied is the distinctive feature of the Mie theory that enables the 
process to be of good use in tropospheric ar.alysis. These ma::ima and 
minima are reduced by increasing the field of view and/or the transmitte 
bandwidth. It is to LID^iR's advantage, then, to retain this complcz 
structure as a data filter. See the Appendix for a sanple plot of C(e<y. 

The complex index of refraction, n, has a dispersion relation lilce 
any electromagnetic medium. See Ta.ble III for an example. 

The less than 1/ change in the index of refraction is t’-nical for 
dielectric media.. The change i' the absorption,y", is even less, otudi 
have indicated that tr.e dispersion in m is balanced by the randomness of 
the aerosol size distribution in such a manner as to almost negate this 
dispersion, 

The angialar dependence is a pre-specified set of conditions, espec¬ 
ially for the case of a LIDiU oachscatter system. This system consists 
of a linearly polarised '.rave undergoing only bachscatter observations. 

As a result of these and similar considerations, one can reduce the 
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number of unlmovr.s drastically by malcing onj.y a fsv simplifying 
assumptions. These ha.ve also been verified by the laboratory analyses 
of many aerosol samples, natural and artificial. 

TABLl III 

Tjiquid water (Q20°G) index dispersion; n = n(l-i')C) 

(A°) n 


4046 .6 

4471.5 

501 5.7 

5460.7 

5292 .6 

6562.3 

7065.2 


1.342724 

1.339423 
1.336363 
1.334466 

1.3329G3 

1.331151 

1.33CC19 


1. The measured values for n raiige from 1.33 to 1.59 '.’ith respect 
to the surrounding air. A large class of aerosols (panticulatcs) may 
be characterised by; 


no = I. 5 4 (i ~ lO. Ooi) 


i'... o, 


Thus, assume that the complex index of refraction is both hnox/n and 
unique. Only the particle’s radius arid number density is allo'/'ed to vary. 

2, Traical non-precipitating aerosols in the troposphere are described 
with their radii bctv;een 0.01^ n and 1 5/«. m. Their number densities can 
be represented by distributions containing only a fe’.; parameters. The 
important model distributions uero presented earlier and the usua:. goal 
of an invers3.on is to best fit one of them. The initial estimates, 
where required, of particulate properties usually includes such models. 

The set of unlmo’-ms reduces to the particle number densities aiul 
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their radii. 3y the use of ore-selected values for r, one only has to 


solve for the distribution functions granting that the radii range is 
properly considered. Therefore the paraaeter of prinary interest is the 
scatterer nuaber density. 

Define the ilie function, K: 

(2S7) 


Uilt; iiJ-c J. Uliu U 

=(2WRYUrr^,^,&) s(^) 


and the particle size distribution function, f(r): 




dhiCrl 

dr 


(2SS) 


The scattering ^xtensity, equation (21 6 ), becoaes: 

K (;<.r) f (O dr (229) 

where x represents the set of user defined parameters (a, X,0, 0). The 
value l(x) is measured, the quantity K(x,r) is computed and the urlaiouTi 
f(r) must be solved for. Tills relation is classified as a Tredholm in¬ 
tegral equation of the first kind. The kernel of the equation is K(o:,r) 
and I(x) represents the transform of f(r) utilizing the stated kernel. 
There exists only a finite set of measurements, l(x), over which the 

'oh. 

inversion may be performed. Let the subscript i represent the i 
measurement of a set of n observations: 

r *"1 

lUO = 3C-itO = 9l = / K (>^i,r) f(r) dr (290) 

' I 

The notation, g^, has been introduced and v/ill be used later in the paper. 
The inversion process is by nature an unstable one and this may easily be 


demonstrated. 


Define: 


sinFr] dr 

Jo 


(291) 


where i= 1, 2, n 

and C is an arbitrary constant. 

(290): 


i: ~ 1, 2, ,,, 


jetting r^ = 0 and r 2 = IT in equation 
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r It 

‘ ^ sin (tr) K (^i ^l')dr (292) 

Since the sin(?r) acts as a weighting function for the kernel on the 
interval (0,11), as ? increases '.-dthout bound, the second integral tends 
to 0. Thus: 

One can anproxinate the transfomed values of f v;ith the function S_, as 
close as possible by talcing ? large enough. Yet the josinYrj function 
nay be nade arbitrarily large and 

?(r) = p(r) + CsLnFr (294) 

will transfom as f(r) with little to no resenblance to f(r) for large 
values of 0, Therefore, one can generate transforas that come as close 
to the desired values, g(x), vdthin any desired accuracy as one may vdsh 
and still have no idea v;hat the accuracy with respect to f(r) nay be. 

The technique used to remove this instability is to require the ad¬ 
dition of a constraining relation on f(r) other than equation (290). 

Tlois nay be done either e:q)licitly or implicitly. The proper relation 
must contain the desired qualities of f(r) requiring a high degree of 
a priori loiouledge. Out of the near-infinite set of possible solutions 
for f(r) generated by the inversion, the constraining relation will act 
as a filter to obtain the unique, or near-unique, solution. 

The presence of measurement error in the g(x^), larger than some 
quadrature error, allows the use of an appropriate quadrature approx¬ 
imation of the form: ^ 

>3 ^ ^ 

for k intervals of the argument y and the w. are the quadrature 

o 

weighting coefficients with increments Ay.;. For equation (290): 
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( 296 ) 


X f (rO A rj 

0-1 ** 

v/here i represents a aeasurenent and j a radii interval. 

Define; for n measurements: 

1, column vector G: i - •••> (297) 

2, coefficient matrix A: a.. = v..K(x.,r.) (29S) 

ij ij ' 

j = 1 ,.,., Ic; Ic radii intervals 

3, column vector II: II. = x(r.) Ar. (299) 

J JO 

Then, equation (296) may he :n:itten as a matrix multiplication: 

G = .III (300) 

The value N. represents the total number of oarticles per unit volume 
J 

’.'.dth radii between r. and r. + ^r., as can be seen via (233): 

- J^. f(r)3r = fCrj) [(D + A rj)- rj] (301) 

and f(r.) is considered constant over the inter/al J. It is further 
J 

assumed that the scattering functions have only small variations with¬ 
in each sub-inter'/al radii range. 

The problem has now been reduced to solving for the k unlcno’.ms in 
II given the (k)x(n) computed values in A and the n measurements in G, 
It must be emphasised that a quadrature approximation v;ill introduce 
serious error in f(r) unless the quadrature v/eighting formula is of an 
accuracy higher than that of the expected measurement error. Let Hq 
be the nu,mber of independent observations made such that n., ^ n. If 
k = Uq, then II may be solved for exactly: 

(require A to be non-singular) 

=r> A''A = i (302) 

I in a matrix equation represents the identity matrix, 

A^J = G => M - A‘'G , (k-Ao') 


where: 

( 300 ): 


(303) 



A 

where 'I represents the interpolated solution vector. However, if k 
n^, then M nay be solved by miniaizing the square of the error (least- 
squares) : 

(require a'^A be non-singular s> a"'"; transpose of A) 

(300): C:r (304) 

this step forns a square matrix which may be inverted, 

M ? ( K nJ (305) 


That this last equation actually minimises the error may be visualised 
by considering an appropriate spatial geometry iraich contains the vector 
G and the vector AH^ where ^ N. See figure 10. The vector is 

seen to be representative of the error vector betv;een G and AH^. The 
minimum error occurs when this error vector has the smallest magnitude, 
thus; 

(AN,) •(<--flN,) =0 (306, 

(the • represents the vector dot product). 

A 

Then for any AH in the range of A, the vector AII^ is still perpendicular 
to it; 

(AA/,).(G-AAi) =o (307) 


’■/riting the dot product as a matrix multiolication requires; 

(flN,)^(G- a 4) = o 

And uy r 0 since it represents any vector in the range, thus; 

(304) 


( 308 ) 

(309) 


In the work that follows, when vnriting the inverse matrix A \ if the 
matrix is not square, the quantity (A.'*’a) ^A"*" is implied. 

The procedure as indicated by equations (303) and (305) are in large 
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n-3pace vector geometry 


G 



Figure 10 


SO 







ooce 


unsuccessful as they do not contain the constraining criterion. 


further problems concern the size of k. 

A 

1, Small k: the vector II does not accurately represent II due 

to the r. set sparsity and/or the crudeness of the quadrature. 

J 

A 

2. Large k; the components of II taice on an oscillatory be¬ 
haviour irrespective of IT and the coefficients in A become particularly 
sensitive to round off error. 

This sensitivity in coefficients is best characterized by the con¬ 
dition number of A (or any matrin) [5] : 

conc/(A) = jflj lf\''l (310) 

■rhere jll - xiora of a matrix 

The norm of a matrix may be defined in any one of three ways depending 
upon the desire of the user to highlight specific as'oects: 

(311) 


1 fli = 

(YICLK 'k 

1 .. 1 


1 1 

1 $ J ^ K *(. r 1 

kcjl 

: niaxirrtum c&LumM - sum 


v;nere 


- magnituae oi tne component 




: maz-timum ro;.-sum (312) 




Iflico = i |nj 

(norriialiy used) |fl]e ~ jui 1^4^ ) *" ■ 3uclidccr. uoru (313) 

Because A"’ must be computed and will contain at least the inaccuracies 
in A, the cond(A) vdll be costly to compute and inaccurate as well. 
IIov;ever, the test is th.c magnitude of the cond(A). If large, the matrix 
A is termed poorly-conditioned and the lil-celihood of a successful inver¬ 
sion decreases rapidly to zero. Oorrections, lilie scaling or partial 
pivoting, must then be applied to A. Ilote, for the identity matrix, 
cond(1) s 1; 

cond (a) > I (314) 

since must exist. Other oronerties of the norm are: 


ol 




lfl[ ^ o 

I fll = O Ip-p /? = O 
1/3/1 I = /3 //? I , K? > O Oind a. numher 
l/l+Bi ^ \R\ +Ib1 
1 B/1 I ^ 161 lAl 


(315) 

( 316 ) 
(317) 
(313) 
(319) 


The ccnd(A) is nost useful in bounding the error expected in the computed 
solution and relating this error to the nagnitude of the residuals. 
Define: 

ill = 3 (320) 

’.■.'here X - exact solution vector 

A 

X - approximate solution *'ector 
3 - vector of hno’-Ti quantities 
A - operator matrix 
The residuals are: 

where R - residual vector 

A 

^ r: _ T 




(321) 

( 322 ) 
(323) 


using ?j bp inverting equation (323) end the norm relations (315)-(319), 
one may incorpora.te the cond(A) as an error measure: 


4 


Rl 


CoMd(A) 161 IXl (324) 

A 

Therefore, the relative error in X, cont 3 ,ined in X, can be as great as 
the relative rosiduaJ. times the cond(A) and as SDcill as this sane re¬ 
sidual divided by cond(A). For large cond(A), the residual p-ields little 

A 

to no information on the accuracy in X. Fquations (320)-(323) may be 
used to define an iterative improvement algorithm (IIA). 
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1 


A A _] 

. compute X by Z - A 


(325; 


(321) 2. compute R by R - E - i-Z 


3. compute 2 by 


_ ,-1. 


( 326 ) 

= R + X (327) 

If the system matrix A has errors in its coefficients, the system 


4. correct X by 


may then become ill-conditioned even if the residuals are small. 

A 

A = A + ?: 

’.rhere A - true coefficients 

H - error in coefficients 
The equation actually being solved is: 

AA 

A r* 'o 


(32S) 


(329) 


(330) 


'7nen, analogous to (324): 

(a) IHt 

Ix' 1A\ 

The relative error in the residual can then be no larger than the rela¬ 
tive error in the coefficients in A times the cond(A). 

Finally, an approximate accuracy check may be obtained from : 

1 £ O'" 






-p 

lO , 


(331) 


then X is probably correct to p digits. 

Equation (300) is assumed to be accurate within a small measurement 
error and equations (303) and (305) are considered insufficient to 
remove the inversion instability. There are several techniques avail¬ 
able to remove thses instabilities of which tliree will be briefly 
covered. 

The first method attempts to smooth the instabilities by one of two 
means. Define: 

a(w) = ((r-AM)\6-ANi) ^ y, (332) 
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where u(i[) - the smoothing measure (a number) 

^1 ’ ^2 ^ 0, numbers 

Mq - a priori information function, vector 
B - matrix describing some desirable smoothing of N 
(G - Pulf (G-AII) - a measure of the accuracy with which N satis¬ 
fies equation (300) 

"ip 

(ri - IIq) (N - ITq) - a measure of the departure of N from the 
expected 

(bn)''* (bn) - a value that measures the departure of N from ideal 
smoothness (BN = O) 

The smoothing technique involves the minimizing of u(N). Some matrix 
differentiation rules will be needed 0-3: 



= zx 


(333) 

AAX 





- A 


(334) 



AZ. 

(335) 

<iX^ 




where X, Y, Z, are vectors and A is a matrix, all of the proper dimensions 
such that the indicated products are defined. 

Condition 1: ^2 - 0; requires a priori knov/ledge, Nq. 


(332): u.Cn)= 


(336) 

where VoC^) " C^'AW} 


(337) 



(338) 

The required condition is: 



= miniraum => U.C’^) 

= 0 . 

M 

(339) 

Utilizing (333)'-('335) ; where 1 is the identity matrix: 


^ = - 2/^’'Cc--an) 


( 340 ) 



(341) 


^4 










Use ( 340 ) and (341) in (339): 

■'■^1 ) (342) 

This is the required solution to equation (300) given the constraint 
function Nq. Note that, by setting = 0, one obtains equation (305). 
Since the multiplier, Vi , represents the difference magnitude bet;;een 
and Nq, one might Oicpect to use the estimated error, (331), as a first 
guess at y, . 

Condition 2: /, =0; requires no pre-conceived solution. 

(332): Uj = Yo M Vt (343) 

v;here - as before (337) 

(344) 

The most common selection for B requires the solution to be smooth 
through the second derivative. In difference equation form: 

Nj/ ^ 


AN 




_ _ M 

•''J ” Arj - rjn-r; - 

^ ^ A Nj ^ Nj*i ~ J_ /Nju ~ Nj i. | Mj+, -/^j 

j ' dr Arj - rjH-rj =Arj ( Arj^i 


(345) 

(346) 

(347) 


Assume equally spaced quadrature points: 

^ CK-I ) 

Tj r JVJI :: Ar 


K - I 


Note thatAr. is independent of j. Thus: 

< =i7*(Mi-2A/iH * Mj,x) 


(348) 

(349) 

( 350 ) 


Since only the coefficients of N. are important, ignore the l/Ar scale 

J 

factor and va*ite: 


N" = (I 


-z 


0 


A 

Mi*i 


85 


(351) 











or U - 3H 

in matrix form where B is a (l:)x(k) matrix defined as; 

/, 


(352) 


1-210 


B = 


0 

0 

0 


1 -2 1 

0 1 -2 
0 0 1 


0 0 0 
0 0 0 
0 0 0 
0 0 0 


V 


0 0 0 0 
0 0 0 0 
0 0 0 0 


1 -2 1 

0 1 -2 
0 0 1 


(353) 


J 


The vector 3II is thus the discrete analog of the second derivative of the 
unlcnoun vector II. Other forms for 3 may be used, notably reference (^20]. 
As before, the constraint shall be to minimise the quantity U2(l0 in an 
effort to achieve the smoothest function II v;ith a minimum of error in Oil: 


UzOO ~ minimum => | ^1 ” ^ 

and Y 2 : 

How, ( 340 ) and (355) in (354): 


(354) 

(355) 

(356) 


This is the required solution to (300) which satisfies the constraints 
outlined above. Again, setting ” 6 yields (305). A recent work 
into this particular solution has yielded considerations for the choice 
of" factor ^('2 is too small, the instabilities arc not 

removed from the final solution. And if chosen too large, the system 
becomes overconstrained and indeoendent of the measurements. The best 
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choice appears tc be ^ sr.d should be so chosen until further 

efforts provide a better choice. 

Both solutions via this technique and the results from the next 
method nay be improved in the iterative fashion as outlined in equations 
(321), (325)-(327); IIA. The algorithm: 

(initial estimate) 1, obtain by (342), (356), or (366) 

2. - G - (357) 

3. (35S) 

4. :r v(°) 4 y(‘'^) (359) 

5. is < some small positive number tT ? 


(m^^' residual) 

(m * error estimate) 
(nev; estimate) 

(test) 


no: m = m + 1, go to step 2. 

''(m+1) 

yes: done, the solution is h to an accuracy on the order 


of 5*, 


The next method involves the computation of some statistical para¬ 
meters. OrJLy the general outline shall be described. See reference [2]} 
for details. Let and be positive definite matrices (all eigen¬ 
values positive) that act as veighting functions, then: 


uhere /7 is (n)x(n) and /Q is (k)x(I:). 


( 360 ) 


Let IIq be the mean vector for N as compiled from previous statistical 
data (model). Define the expected value operator E as: 

l,rna 3 (X-,) -PifkO ( 351 ) 

where f(x) is the probability density function for g(x). Let the vector 
£ have zero mean and measure the error in G. Define the covariance 
matrices thus: 
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5. = E&E-3 


(363) 


Assuning a nornal distribution for 6 and 'I, then the Gauss-Marico’ 


theoren yields a statistically optiaun estimate II by minimizing v(Il) 

s 


for: 


r, - 5, 


-I 


(364) 

(365) 

( 366 ) 

(367) 

A _1 

The resulting covariance matri:: for (II - II ) is found to be X . Thus. 

s 


then: 
where 


G = 

D, = x''a"5;'(c>-aMoVN, 


• 3 - /N r, W^g_ 

X = S'* + A 


the expected mean square error may be ’.vritten as: 

= trace X ' 


(363) 


I-Iithout solving equation (300), one can thus obtain an accuracy estimate 

A 

for il under the stated statistical conditions. One can also annly the 

IIA, (357) - (359), to enhance the estimate II . 

s 

The last teclinique to be considered is an iterative scheme somev/hat 
analogous to the iterative Improvement algorithm. The method used is 

r 1 ^ ’ 

the Land^;eher iteration*jj^y. La general, an iterative giiessj 'L ' is 
made ba. 3 ed on sone a. priori IcnoT/ledge and then inproved successively 
to form a sequence ’..'hose member, I. , transforms to the measured 


(m ] 

' for smaj.l m. Tills 


G as m tends to infinite. The techniques' s’uccess is due to the auto- 
ma.tic constraJnin.g of the oscillatory nature of II* 
ahlcws an approximation to C- early enough to avoid the oscillations if 
the secuence is terminated when the residual error is cn the order of th 


error in G, The Landvreber algoritiun uses; 


(369) 


’,/nere 


hlC) 


reflects the eic'ected solution. If the coefficient matrix 


S8 












l:as "j-l _?osi'biv9 cor._-onent3, the conver^^ence i/ill 
finii's e. diagozial i^atrirc D: ^ 

" Cl-, j=I."*, 

arid reva’itir." the estimate as: 

Cm) A Ci>''0 «T/'/' n ifi 

=N^ + tiA’^CG-flWo ) 


be ir."ro'/8d b’’’ ds- 




(371) 


This conpletes this short survey of inversion techniques. A coanon 
feature for ell of thee is their use of a priori Imouledge to constrain 
the selection of probable solutions in order to enhance the choice of 
the proper solution. Indeed, this criterion nust be net by any tech¬ 
nique the ezroerinenter chooses for inversion, if he is to have anp' 
confidence in the results. 


' r Tto " 


This paper has endeavored to present the requirenents for a device 
designed to nap the tropospheric aerosol properties on a real-tino basi 
The specific ncdel equations ’..-ere derived a-fter a theoretical devclopne 
■;a.3 nerforned in such detail the.t the user becane familiar ’.;ith bhe 
limitations of the theory, A candidate for such a device that uouj.d 
utilise advantageously the specific thscry characteristics is the laner 
radaa-, or lIDfl’,. The LID.'l systen allo\<s for the direct selection of 
the hey pananaters of the theory and yields as output the required in- 


:uts for the data inversion to bevin. A feu Inversion teclmlaues were 


presented, but, due to the lou’ge ru.V’.ilaAl.e entent of such net’nods, the 
subject was by no means exhausted. The aerosol LIDPJ. must be capable 
of polanisation measurements and fast wavelength scanniing to increase 
the num'oer of indeoendent observations made. 3uch system ca.nability 











\/ould decrease tlie rceuired arour.t of a "oriori laior.'ied-e the eineri- 


nsntcr r.ust decide on. 

foth nresont and future developnents 
lasers argue strongly for the design, co 
aerosol analysis system. 


in fast processors 
nstriaction, and use 


and versat 
of a LIDif?. 


\ 


90 















































































































;u-P 3IDD:: CCI!?UTATI0:L\L jX)R:-rJLx\2- i\in3 :\LQCRITH'I 


This section presents sone further equations to facilitate the conpu- 
tation of the required fuctions fron the Mie theory. 

Bachscatter: 0 = 1c0° 


Tr„ (i8oy = (- 1 )''"' f 

TnOso”) = -TT^ (iW”) 

5,(l9£)’} = ^ , C-0"(n )( b„-an) 

3i(l80’) = -S, (l80°) 


(372) 

(373) 

(374) 

(375) 


backs c a’ 


1 I 2 

tter gain: (>= ^ , C*0 Cn+l)(b^-o.n')| ( 376 ) 


scattering functions: 


recurrence relations: 

7^ (cos e) = C05 © Tin (‘^ 0 -^ ' sin^eiTn (cos e) (377) 

■rr„ (C 05 e) = (cos 9)(^) ir„., (cos o) -(^) 17., (cos e) ,,,,, 

■n'nf6£)5o) = f2'’-')rr„.,(cose)+ft-42(co5 e) ( 379 ) 

1 ("o (C05 ©) = o (3S0) 

TT, (cos 0) = I (331) 

'n;'(co 5 e)-o (322) 

1T'/(C0 S©)= O (3S3) 


Scattering coefficients ; The equations for a^ and may be vn^itten in 
terms of the logarithmic decrement function,(Z). The equations 
below have been so x-nritten and then incorporated into an algoritrim that 
was used to compute the scattering and extinction efficiencies for a 
particular case of m = 1.54(l-i0«C0l), These values vfere then plotted 
on the enclosed graphs to picture the oscillatory behaviour of the 
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scattering efficiency, Q , for <=^ fron 0«1 to 17.0 in 0.1 increments. 
The data points were connected by straight line segments to complete the 
visualisation. It must be realised that further structure may appear 
under finer resolution. The routine was programmed for the heulett- 
Pachard model 410 machine vrith the program listing accompanying the 


figures. 



(384) 

m : ^optC'-■<•><') = 

(385) 

/3 - mcx = /3>c•^ ^/3y 

(3S6) 

llote in the follouiig that orJLy ^ requires 

a complex argument. 

(iSo) - -ij |_^C4n+0 

(187) QscaC'*^) - iy + 1 1 



(387) 

, y. m 1 C/s)_ 1 

'='> -■ ter-;)- -”7? ('3) J 

(383) 

yj (or j = sio ^ 

(389) 


(390) 


(391) 

Xo C°<) = C 05 Coc) 

(392) 

X, (^) = 2^^ + % 

(393) 

x„,M= c^x„C-)-x„.,te 

(394) 

(114) Jn(^)= tW+iXnC®^') 



(395) 

x;;te= X„(■-;) 

(396) 


(397) 

i'V) - 3n C”<2 

(398) 

(399) 
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/ 


(400) 


tar? Z^* 


-Al.^crithn 


=> Qsca(f“) = ^ %' (^) 

%t(=<)= (2n ♦ i) f?c fa, C«<) <■ t„ C»<) j 

=> QextW z: ^'?2.Z«) 

Begin: A, 1, input: opt -i ^ 

2. coapute: rrtx ^ y 

B. 1. input: (S^cT (required accuracy) 

2. compute; 5 y 

3. assign; n^j C.^) - %z(^) - O 

4. compute; 

r.„c^) 5 X,c»<) 

f/3) 

C. 1. compute: ("oc) ^ ^6><) 

C”-) 


?0) 


2. compute; CL^, ^ 

3. conpute: ? = ^, f-) + (2 « ►<) f |a„ T^)) t / C^) P} 


= gz («) t-(zr,ti ) Re fa,C»‘) + b, 0*) ] 


4. test: is 


c*=<) 

yes: go to D.1, 
no: go to C.5. 


5. test: is K(,"V°<)-<^y^)) ^ ^ ? 


^ C“) 
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yes: go to D.l. 
no; go to S.1. 


D. 1 


1. 

assign; 

11 






2. 

conpute; 


X. c 

3. 

assign; 

V'„., (“) 

= 'f'r, ('=') 



/„ = 




y^n-l (<^) 




Xn M = 

('=^) 




= 7'.'' 



n = n 1 1 


4. 

go to G.' 

1. 


1. 

conpute: 

*^sca 



%xt C'^) " 


2. 

output: 

n(nunber of 

steps) 


Q (^) 

sea 

*^ext (oi') 
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attcrlng efficiency versus size parameter for sphere 
m = 1.54( 1-iO.OOl) 


O 

CO 
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Size parameter 
























cattering efficiency versus size paraxnetor for slipere 

m = 1.54(1-10.001) 
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